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1.  Introduction 

This  is  the  final  report  for  Contract  F19628-90-C-0097,  “Radiation  Belt  Dynamic 
and  Quasi-Static  Modeling  Based  on  CRRES  Data.”  The  period  of  performance  covered 
by  this  final  report  is  1  May  1990  through  30  June  1994.  The  research  and  modeling 
efforts  conducted  under  the  contract  have  achieved  the  essential  elements  of  the  original 
objectives.  However,  due  to  the  unfortunately  premature  demise  of  the  CRRES 
spacecraft,  we  were  able  to  coUect  only  a  minimal  data  base  to  test  quantitative  elements 
of  the  dynamic  modeling  effort.  Nevertheless,  for  short-tenn  responses  of  the  outer 
electron  belt  to  geomagnetic  disturbances,  we  were  able  in  each  dynamical  event  to 
delineate  quantitatively  the  extent  to  which  the  relativistic-electron  response  was 
“adiabatic”  and  the  extent  to  which  it  was  “diffusive”  (usually  corresponding  to  a  loss  of 
particles  from  the  magnetosphere).  For  the  intermediate-term  response,  we  were  able  (by 
using  a  SCATHA  data  base  organized  with  respect  to  pitch  angle,  L  value,  and  energy)  to 
constmct  and  verify  a  coir^rehensive  model  of  simultaneous  difiusion  of  radiation-belt 
particles  in  pitch  angle  and  L  value,  taking  account  of  the  energy  change  that  inherently 
acconpanies  radial  difriisian.  In  addition  to  data-based  modeling  of  the  dynamic  radiation 
belt  with  CRRES  and  SCATHA  data,  we  were  able  to  set  up  single  particle  trajectory 
tracing  codes  as  ancillary  tools  to  understand  the  dynamic  behavior  of  energetic  radiation 
belt  particles  “contained”  in  the  “leaky”  magneto^here.  Although  we  were  not  able  to 
exercise  these  tools  extensively  on  the  CRRES  data,  due  to  curtailment  of  fimding 
subsequent  to  satellite  failure,  tfiey  nevertheless  prove  to  be  valuable  for  modeling  the 
outer  magnetosphere.  Technical  summaries  and/or  literature  citations  to  each  of  the  above 
achievements  will  be  given  in  this  report. 

In  the  rest  of  the  Introduction  to  this  final  repon,  we  faighligfat  the  history  of  the 
project  and  its  original  objectives.  Summaries  of  substantial  achievements  imder  the 
contraa  are  given  as  separate  technical  sections  of  this  final  report.  Finally,  a  publication 
list  containing  six  items  (three  articles  in  refereed  scientific  journals,  two  papers  in  refereed 
monographs,  and  one  journal  manuscript  recently  submitted,  aU  directly  attributable  to  this 
contract)  is  also  furnished. 

Existing  characterizations  of  the  natural  radiation  belts  lack  sufficient  accuracy  to 
provide  for  reliable  trade-off  studies  between  space-system  survivability  and  other 
spacecraft  systems  insofar  as  short  and  intermediate  term  radiation  belt  responses  are 
concerned.  Lx)ng-term  fluences  of  penetrating  electrons  and  protons  obtained  by 
integrating  the  current  NASA  flux  environments  in  typical  low-Earth  and  geosynchronous 


Dibits  are  considexed  accurate  only  to  within  a  factor  of  about  two  (up  or  down). 
Fluences  in  oibits  throu^  other  regions  of  the  magnetosphere  are  even  more  unceitain, 
and  individual  flux  listings  may  be  in  error  by  more  than  an  order  of  magnitude.  This 
major  inadequacy  in  environment  specification  often  expresses  itself  in  the  form  of  severe 
shielding-weight  penalties  imposed  (xi  suivivable  system  payloads,  or  as  unforeseen  human 
risks.  Large  and  conplex  ^ace  systems  can  ill  afford  these  lisks,  which  are  paiticulaily 
severe  when  humans  are  involved.  The  weight  penalties  for  flion-duration  payloads  may 
be  even  worse;  lack  of  predictive  aq)ability  means  that  they  must  qualify  against  the  most 
severe  environments,  which  may  not  even  be  relevant  over  a  mission  of  short  duration. 

This  project  was  amceived  in  response  to  the  dual  diallenge  of  constmcting 
dynamic  models  of  the  radiation  belts  and  of  supporting  the  development  of  quasi-static 
models.  This  effort  is  specifically  based  on  the  coii^rehensive  database  obtained  from  the 
cross-calibrated  conplement  of  particle  radiation  instruments  provided  by  Lockheed  and 
by  AFGL  on  the  CRRES  and  S<IATHA  missions. 

Radiation  belts  are  presently  characterized  in  terms  of  montages  of  satellite  data 
san:q>led  over  various  spatial  regions  regardless  of  time,  with  uneven  statistical  distribution 
and  with  only  rough  categorization  in  terms  of  geomagnetic  activity.  We  describe  the 
compilation  of  such  data  moitages  as  static  modeling. 

The  next  level  of  radiadon-belt  characterization  is  achieved  by  recognizing  the 
importance  of  dyiuonic  changes  of  particle  fluxes  in  response  to  geoiiuignetic  activity  or 
other  stimuli.  This  level  seeks  to  model  tenoral  variations  of  the  particle  distribution 
function/  according  to  a  set  of  a  priori  physical  laws,  as  ^lied  to  an  assumed  model  for 
the  geomagnetic  field.  A  satellite  database  is  used  to  define  flie  relevant  geophysical 
parameters  on  the  one  hand,  and  to  verify  the  adequacy  of  the  model  assunq)tions  on  the 
other.  We  describe  the  physical  characterization  of  the  evolving  phase-space  distribution 
function,  in  terms  of  geomagnetic  parameters,  as  dynamic  modeling.  It  should  be  kept  in 
mind  that  dynamic  modeling  seeks  to  achieve  a  measure  of  space-time  pnedictive  C29)ability 
through  identification  of  physical  pnocesses  and  geomagnetic  specifications;  it  is  not 
necessarily  relevant  to  the  reduction  of  uncertainty  in  static  models,  as  hose  uncertainties 
are  not  necessarily  of  physical  origin. 

Even  at  the  level  of  impioving  statistical  data  montages  (static  nnodels)  to  include 
categorization  geomagnetic  states  (quasi-static  modeling),  the  supjprort  of  dynamic 
modeling  cannot  be  dispiensed  with.  For  a  single-satellite  constellation,  he  geomagnetic 
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state  is  likely  to  have  changed  before  a  complete  spatial  saiii(ding  of  data  can  be  collected. 
This  is  paiticulaily  troublesome  for  very  quiet  or  very  active  states  of  the  magnetosphere, 
as  the  IMF  is  very  variable  during  such  periods.  Incomplete  spatial  sampling  over  a 
specified  temporal  state  can  totally  undermine  the  statistical  basis  of  data  montages.  Since 
dynamic  modeling,  after  verification,  provides  a  means  of  resolving  spatial,  temporal,  and 
magnetic  effects,  it  provides  (in  our  view)  the  rally  means  of  supporting  the  upgrade  of 
statistical  data  montages  into  quasi-static  models  keyed  to  geomagnetic  parameters. 
Dynamic  modeling  is  the  major  objective  of  our  effort  and  is  recognized  to  be  a  major 
unsolved  scientific  problem. 

The  primary  objective  of  the  contract  was  to  construct  a  dynamic  model  of  the 
radiation  belts,  a  model  dependent  on  phase-space  variables  and  time,  based  mainly  on 
measurements  of  the  ONR-307'3  SEP  instrument  onboard  (ZRRES  and  on  constraints  set 
by  the  physics  of  particle  transpon  in  the  magnetosphere.  Other  relevant  CURES  data 
have  been  included  subject  to  availability.  Another  primary  objective  was  to  eiqiloit 
dynamic  modeling  techniques  to  support  the  construction  and  inqnovement  of  quasi-static 
radiation-belt  models.  The  specific  objectives,  each  of  Miicfa  is  referenced  to  a  proposed 
task,  are  listed  below: 

1.  Reduce,  refine  and  organize  SEP  database  for  (tynamic  modeling.  Acquire  other 
(TRRES  data  subject  to  availability  and  suitability. 

2.  Test  and  improve  existing  simultaneous  bimodal  diffiisirai  model  with  the  model 
database,  extending  its  current  applicability  in  the  outer  electron  belt  toward  the  inner 
electron  and  proton  belts.  The  energy  range  of  concern  is  >  200  keV  for  electrons  and  >  1 
MeV  for  protons.  The  resulting  dynamic  model  will  beconoe  the  AFGL/Lockheed 
Dynamic  Radiation-Belt  Model. 

3.  Model  physical  j^enomena  that  need  to  be  incorporated  into  the  AFGL/Locldieed 
dynamic  nx>del.  These  include  shell  splitting,  effects  of  encounters  with  the 
magnetopause,  charge  exchange  with  the  atmosphere,  and  wave-particle  interactions. 

4.  By  examination  of  the  database  and  model  results  for  a  sample  of  no  less  than  20 
magnetic  storm  episodes,  constmct  period-averaged  snqrshot  models  of  the  storm-time 
radiation  belts  for  the  purpose  of  clarifying  the  stimulated  response  of  the  radiation  belts. 

5.  Dynamic  modeling  techniques  will  be  exploited  to  support  the  upgrade  of  static 
models  into  quasi-static  rtKxlels.  A  corollary  benefit  of  dynamic  modeling  is  that  it 
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provides  a  method  for  projection  of  data  obtained  on  a  sateHiic  trajeaoiy  to  off-trajectoiy 
space-time  points,  beyond  the  sin^  mapping  of  fluxes  along  a  magnetic  flux  mbe.  Ihis 
technique  can  be  tqiplied  as  a  screen  to  reduce  the  lisk  of  combining  data  from  diverse 
magneto^heiic  states  in  the  consoucdon  of  quasi-static  models. 

6.  Investigate  the  feasibility  of  invariant-space  difiuskn  as  a  second-generation 
develc^jment  of  radiatian-belt  dynamic  modeling. 

As  will  be  shown  in  the  teduiical  sections  below,  we  have  achieved  and  published 
most  of  the  scientifically  relevam  elements  of  these  objectives  with  various  degrees  of 
depth,  despite  the  downscqiing  of  the  effon  after  the  eafiy  demise  of  CRRES.  Evidence 
of  our  achievements  is  furnished  in  the  following  technical  sections  (Sections  2-5)  and  in 
the  list  (Section  7)  of  publications  that  have  resulted  from  this  effort. 

2.  Relativistic-Electron  Response  to 

Intensities  of  geomagnetically  mqjped  outer-zone  radiation-belt  electrons  typically 
decrease  (often  by  orders  of  magnitude)  as  the  main  phase  of  a  magnetic  storm  develops, 
then  increase  at  least  bade  to  pre-storm  levels  (and  often  higher)  during  recovery  j^iase. 
The  “deanest”  events  are  those  in  which  a  single  decrease  in  (from  values  near  zero 
toward  values  around  -100  nT)  is  followed  by  a  smoodi  recovery  without  further 
fluctuation.  An  inqxrrtant  part  of  the  work  performed  under  die  present  contract  was  a 
systematic  study  of  the  typical  r^ladvistic-election  response  to  such  changes  in  D^.  This 
study  invoked  a  CRRES  data  base  induding  about  15  of  such  “clean”  events  (as  well  as 
some  less  “clean”  events)  that  occurred  during  the  lifetime  of  CRRES.  The  purpose  of  the 
smdy  was  to  detetmitK  the  extent  to  vfrich  the  typical  outer-belt  electron  response  to 
is  essentially  adiabatic  (consistent  with  consovation  of  the  diree  adiabatic  invariants 
during  a  global  inflation  and  subsequent  relaxation  of  drift  didls)  and  to  what  (residual) 
extent  the  typical  sttnmtinK  decrease  in  energetic-electron  intensities  must  involve  a 
genuine  loss  of  mqiped  partides. 

The  study  invdeed  a  novel  method  of  analydng  energetic-electrmi  and  magnetic- 
fidd  data  from  CRRES,  so  as  to  reveal  the  spatial  structure  of  the  magnetic-field 
perturbation  AB  needed  in  tnder  to  account  quantitatively  for  the  main-phase  decrease  in 
particle  flux  as  Dg^  decreased  toward  more  negative  values,  hi  several  cases  the  spatial 
structure  thus  required  bore  a  dose  resemblance  to  the  radial  variation  of  AB  anticipated 
from  standard  ring-current  models  ,  although  not  to  the  radial  variation  of  the  AB  actually 
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observed  in  those  cases  (see  tqipended  preprint  for  details).  A  more  traditional  analysis  of 
the  stomitime  energetic-electron  response  to  changes  in  based  on  the  adiabatic 
transformation  of  radial  profiles  and  energy  spectra  observed  between  storms,  suggests 
that  adiabatic  response  can  account  for  anyvt^iere  from  zero  to  100%  (but  typically  less 
than  half)  of  the  observed  stoimtime  logaridunic  decrease  in  electron  flux,  but  that  scatter 
among  the  data  points  makes  it  difficult  to  make  »  more  quantitative  estimate  than  diis. 

A  less  soohisticated  analysis,  based  solely  on  a/tiahatir  transformation  of  observed 
pre-stonn  energy  spectra  in  conjimction  with  a  model  (see  appended  preprint)  for  the  AB 
produced  by  the  ring  current,  suggests  that  adiabatic  energy  loss  might  reduce  the  trapped 
outer-zone  electron  flux  at  L  *  4.5-6.5  by  a  factor  ~  2  for  =  -35  nT  or  1^  a  factor  - 
10-15  for  =  -135  nT.  Reduction  factors  actually  observed  in  the  CRRES  data  base 
ranged  from  about  10-50  for  Dg^  =  -35  nT  to  about  100-400  for  =  -135  nT.  This 
comparison  suggests  that  adiabatic  energy  loss  alone  is  usually  inadequate  to  account  for 
more  than  30-50%  of  the  typical  stoimtime  logaiiflimic  decrease  in  energetic-electron 
fluxes  at  L  ^  5,  even  if  the  upsetting  flux  increase  associated  with  outward  motion  of  the 
radial  profile  is  neglected.  The  typical  flux  drop-out  thus  appears  to  represent,  at  least  to 
some  extent,  a  genuine  loss  of  trapped  energetic  elections  frixn  die  magnetosphere. 

Such  an  inference  is  not  entirely  surprising,  as  relativistic  electron  precpitahon 
(REP)  has  long  been  known  [Forbush  et  al.,  1961]  and  understood  [Thome  and  Kennel, 
1971]  to  accompany  the  main  |diases  of  geomagnetic  storms.  Indeed,  a  probabilistic 
model  for  the  outer-belt  electron  intensity  [Chiu  et  al.,  1980,  pp.  137-140]  is  based  on  the 
stormtime  occurrence  of  this  phenomenon.  However,  the  present  work  constitute^  a 
careful  study  of  the  relative  contributions  of  adiabatic  response  and  electron  piecpitation 
to  the  observed  stoimtime  modulation  of  relativistic  electron  intensities. 

A  manuscript  describing  these  findings  (Short-Term  Responses  of  Outer-Belt 
Relativistic  Electrons  to  Variations)  has  been  submitted  to  the  Journal  of  Geophysical 
Research.  A  copy  of  this  manuscrpt  [Schulz  et  al.,  1994]  is  attached  as  an  appendix  to 
fliis  final  report. 

3.0  Particle-Trigectory  Computations 

The  purpose  of  paiude-trajectoiy  tracing  in  work  performed  under  this  contract  is 
to  elucidate  certain  transport  processes  for  which  diffusicHi  coefficients  are  difficult  to 
calculate  in  the  usual  way.  The  usual  perturbation  theory  of  charged-particle  transport  is 
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based  on  the  tnodificadcm  of  relatively  siin{de  adiabatic  trajectories  by  small  perturbations, 
such  that  violation  of  orre  or  more  of  the  adiabatic  invariants  can  be  described  a 
diffusion  coefficient  proportional  to  the  square  of  the  amplitude  of  the  perturbation.  The 
present  transport  investigation  involves  a  violation  of  one  or  more  adiahanr  invariants  in 
the  course  of  particle  motion  through  the  urq)erturbed  magnetosphere.  The  process  of 
interest  here  entails  a  violation  of  the  secoixl  adiabatic  invariant  at  the  dayside  bifurcation 
of  the  magnetic-equatorial  surface  (locus  of  minima  in  B  along  field  lines).  The  bifurcation 
affects  guiding-center  trajectories  that  jq^proach  within  about  2  /{g  of  the  dayside 
magnetopause. 

Expected  consequences  of  this  equatorial  bifurcation  for  the  adiahatic  motion  of 
charged  particles  had  been  outlined  by  Shabansky  [1971].  In  particular,  the  second 
invariant  of  a  particle  approaching  the  bifurcation  from  the  unbifurcated  side  has  been 
expected  to  partition  itself  between  two  sub-invariants  corresponding  to  smaller-amplitude 
oscillations  about  the  northern  and  southern  branches,  respectively,  of  the  minimum-B 
surface.  (This  works  out  to  a  half-and-half  partition  only  when  die  magnetic  d^le 
moment  is  perpendicular  to  the  solar-wind  velocity.  (Mierwise  the  partition  is 
asymmetric.)  Conversely,  the  two  sub-invariants  have  been  expected  to  reven  to  the 
original  second  invariant  upon  traversal  of  die  bifurcation  from  the  bifurcated  side. 

Particle  trajectory  computations  undertaken  in  the  present  smdy  do  not  entirely 
confirm  the  anticpated  quasi-conservation  of  the  second  invariant  Violation  of  the 
second  invariant  (and  also  die  diird  invariant)  is  especially  perceptible  for  particles  that 
mirror  near  the  magnetic  equator. 

Simulations  performed  in  connection  with  this  contract  were  made  in  the  “source- 
surface  model”  of  the  magnetosphere  {Schulz  and  McNab,  1987].  This  is  a  generalization 
of  the  solar-coronal  magnetic-field  model  of  that  name  [Schatten  et  al.,  1969].  The 
magnetospheric  source-surface  model  is  a  “prescribed-magnetopause”  model  like  those  of 
Voigt  [1981]  and  Stern  [1985].  It  admits  an  almost  arbitrarily  specified  magnetopause 
stupe.  The  source-surface  model  provides  for  a  geomagnetic  tail  in  \«^ch  the  magnetic 
field  lines  are  constructed  a  priori  by  analytic  geometry.  In  its  present  form  the  model 
accommodates  an  arbitrary  ang^e  \|r  between  the  solar  wirxl  velocity  (a)  and  the  planetary 
magnetic  moment  (p,).  The  case  Y  =  90®  was  treated  by  Schulz  and  McNab  [1987].  Any 
realistic  description  of  diurnal  and/or  seastmal  variation  in  the  Earth’s  magnetosi^te 
requires  an  extension  to  arbitrary  values  of  \|r,  and  this  turns  out  to  have  been 
straightforward.  For  example,  the  source-surface  model  yields  the  magnetic-field 
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configuration  illustrated  in  Figure  1  [Schulz  and  McNab,  1994]  for  y  =  60°.  The  distance 
b  from  the  point  dipole  to  the  nose  of  die  magnetqianse  is  taken  as  10  for  purposes  of 
illustration,  but  this  is  likewise  an  adjustable  parameter  of  the  model.  The  left-hand  panel 
in  Figure  1  shows  field  lines  that  emanate  from  a  d^xile-centered  sphere  of  radius  1  at 
5°  intervals  of  magnetic  latitude  in  the  noon-midnight  meridional  {dane.  The  source- 
surface  model  yields  (as  ou^ut)  an  unambiguously  defined  and  realistically  warped  neutral 
sheet  whenever  the  angle  y  between  the  the  Earth’s  dqiole  moment  p  and  the  solar-wind 
velocity  u  is  ^  90°.  The  right-hand  panel  of  Figure  1  shows  the  intersection  of  this  neutral 
sheet  (solid  curve)  and  of  open  field  lines  (chosen  so  as  to  emanate  from  the  above- 
described  dipole-centered  sphere  at  5°  intervals  of  magnetic  latitude  and  15°  intervals  of 
magnetic  local  time)  with  a  distant  downstream  cross-section  of  the  Earth’s  magnetic  tail 
for  y  =  60°. 

An  especially  relevant  exan^rle  of  trajectory-tracing  (in  this  case  for  y  =  90°)  is 
shown  in  Figure  2,  where  an  “equatorially  minoring”  proton  is  started  on  the  night  side 
(around  0200  MLT)  and  allowed  to  move  in  accordance  with  the  Lx>rentz  force  equation. 
The  motion  resembles  a  clockwise  guiding-c^iter  drift  until  about  1600  MLT,  where  the 
particle  encounters  the  well-known  bifiirc^on  [Shabansky,  1971]  of  the  magnetic 
equatorial  surface  (locus  of  minima  in  s  IBI  along  field  lines).  The  proton  had 
meanvdiile  been  bouncing  with  small  amplitude  and  evidendy  its  ’’guiding  center”  had  a 
slightly  positive  Vn  at  the  time  it  encountered  the  bifurcation,  llius,  it  chose  to  follow  the 
northern  branch  of  the  equatorial  surface  (dense  trace,  right-hand  panel).  Shabansky 
[1971]  showed  that  the  second  adiabatic  invatiant  J  is  ideally  paititioned  at  the  bifurcation 
site  (and  thus  cut  in  half  if  the  B  field  is,  as  here,  symmetiic  between  north  and  south). 
However,  the  flat  and  curved  portions  of  the  mininaim-B  surface  intersea  sharply  (at  right 
angles  in  the  symmetrical  case)  at  the  bifurcation  site  [e.g.,  Whipple,  1979],  and  the 
second  adiaharic  invariant  qrpears  to  be  violated  as  a  result  The  geometry  of  the 
minimum-S  surface  requires  an  abrupt  change  in  direction  for  the  guiding-center  drift 
motitHi  for  particles  that  mirror  near  the  magnetic  equator.  Evidently  the  required  change 
in  direction  is  just  too  abrupt  to  maintain  7  =  0.  The  bounce  amjditude  in  Figure  2  seems 
to  have  increased  as  the  proton  started  to  fcdlow  the  nortiiem  branch  of  the  minimum-^ 
surface  around  1600  MLT  and  again  as  it  exited  the  nortiiem  branch  around  0800  MLT. 
The  particle  thereafter  bounced  visibly  between  northern  and  southern  mirror  points,  still 
executing  a  noi-chaotic  motion,  until  it  almiptly  escaped  into  die  tail.  (If  the  first 
invariant  is  conserved,  then  violation  of  the  second  invariant  typically  leads  to  violation  of 
die  third  invariant  also,  in  this  case  by  enough  for  the  third  invariant  to  lose  its  meaning.) 


7 


8 


9 


In  other  examples  of  simulated  panicle  moticn  [e.g.,  Chiu  et  al.,  1991,  Figs.  33- 
35,  rq)roduced  here  as  Figs.  3—5],  such  that  the  representative  particle  had  initially 
mirrored  quite  far  off  the  equator,  no  such  violation  of  the  second  or  diird  adiffhatic 
invariant  had  been  evident 
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Figure  3.  Projected  motion  of  1-MeV  proton  started  with  local  pitch  90.9°  at  (X,  F,  Z)  = 
(7,  3, 0)  in  conq)ressed  version  of  paraboloidal  model  magneto^here.  This  is  a  projection 
onto  the  equatorial  plane.  Starting  point  is  in  equatorial  plane  but  wdl  away  horn  i^itber 
branch  of  bifurcated  mmimum-fi  siMace. 
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motion  of  1-MeV  proton  started  with  local  pitch  90.9°  at  (X,  Y,Z)  = 
(  .  >  ^  compressed  version  of  paraboloidal  model  magnetosphere.  This  is  a  projection 

onto  the  noon-midnight  meridional  plane.  Bounce  motiiHi  is  visibly  asymmetric  between 
north  and  south  because  minimum-B  surface  is  bifurcaied  near  dayside  magnetopause. 
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Figure  5.  Projected  moaon  of  1-MeV  proton  staned  with  local  pitch  90.9®  at  (X,  y,Z)- 
(7, 3, 0)  in  compressed  version  of  paraboloidal  model  magnetosphere.  This  is  a  projection 
Goddng  from  the  Sun)  onto  the  dawn-dusk  meridimal  plane.  Bounce  motion  is  again 
visibly  asymmetric  between  noith  and  south. 
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4.0  Diffusion  Studies 


Major  efforts  were  devoted  during  the  eaily  phases  of  diis  contract  to  an 
examination  of  the  behavior  of  trapped  hi^-eneigy  magnetospheric  particles  under 
violation  of  one  or  more  of  the  adiabatic  invariants.  Such  particles  were  investigated 
collectively  through  CRRES  (and  other  satellite)  data  imeipretation  designed  to  identify 
the  radial-diffiision  and  pitch-angle  diffusion  coefficients  simultaneously.  Such  particles 
were  also  investigated  individually  through  numerical  siinulations  of  dieir  motion  in  the 
outer  magnetosphere,  where  the  conditiems  for  magnetic  contaiiunent  are  not  well 
established.  Some  of  the  results  of  this  work  have  been  puUished.  Other  results  were 
summarized  in  the  interim  report  (USAF  PL-TR-92-2160).  Rather  than  reproduc  this 

material  here,  we  cite  the  relevant  publications  and  provide  only  brief  summane  the 

diffusion  studies  here. 

4. 1  Simultaneous  Radial  and  Pitch- Angle  Diffusion  in  die  Outer  Electron  Radiation  Belt 

A  solution  of  the  bimodal  (radial  and  pitch-angle)  diffusion  equation  for  the 
radiation  belts  was  developed  with  special  regard  for  the  requirements  that  pertain  to  the 
analysis  of  radiation-belt  data  acquired  from  satellites.  In  die  p^r  bearing  the  above  tide 
[Chiu  et  ai,  1988]  we  used  diis  solution  to  test  the  bimodal-difiusion  theory  of  the  outer- 
belt  electron  distribution  by  confronting  the  solution  with  satellite  data.  Satellite 
observations,  with  coverage  inhetendy  limited  in  L  and  dine  (r),  sddom  provide  sufficient 
coverage  of  phase  space  to  describe  relaxadon  of  the  entire  radiation  belt  to  equilibiium. 
Indeed,  the  radiation  belt  as  a  whole  may  never  attain  equflihiium.  Moreover,  since  time 
intervals  for  i^ch  continuous  data  coverage  is  available  are  often  comparable  in  length  to 
geomagnedc  disturbances,  it  would  be  inapprt^iriate  to  interpret  such  data  (from  a 
window  limited  in  L  and  t)  in  terms  of  a  diffusion  theory  based  on  semi-infinite  temporal 
responses  to  inqiulsive  disturbances.  The  observadonal  limitadons  diemselves  ini|dy  that 
qipropriate  solutions  for  the  interpretation  of  satellite  data  are  general  solutions  for  a 
finite-volume  boundary-value  problem  in  bimodal  difiuskn  theory.  In  the  cited  piper 
[Chiu  et  al.,  1988b]  we  qrplied  such  a  solution  to  the  analysis  of  radiation-belt  electron 
data  acquired  by  the  SCATHA  satellite  during  moderate  geomagnetic  activity  and  thoeby 
fouixl  that,  because  of  its  generality,  this  ^rproach  is  a  promising  candidate  to  be  the  basis 
for  a  new  ^rproach  to  the  dynamic  modeling  of  radiation  belts. 
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4.2  Toward  Dynamic  Modeling  of  the  Outer  Electron  Radiation  Belt 


In  the  paper  bearing  diis  title  [Chiu  et  al.,  1990a]  we  showed  that  the  general 
solution  of  the  simultaneously  bimodal  (radial  and  pitch-angle)  difhision  equation  can 
represent  the  SCATHA  outer-belt  electron  distribution  to  a  hi^  degree  of  accuracy  over 
time  intervals  up  to  about  10  hours.  Successhil  rq}resentation  of  the  data  was  found  to 
require  diffusion  coefficients  with  parameters  having  numerical  values  entirely  consistent 
with  the  range  obtained  in  previous  observational  and  theoretical  studies.  Such  a 
representation  of  satellite  data  satisfies  the  basic  requirements  of  a  dynamical  model  for 
the  outer  electron  belt  at  L  <  7  under  both  quiet  and  disturbed  geomagnetic  conditions. 
For  L  ^  7  the  representation  encounters  difficulties  on  account  of  “butterfly”  pitch-angle 
distributions,  which  mr^  indicate  severe  ^lell  splitting  effects  or  even  energetic-electron 
encounters  with  the  magnetopause. 

4.3  Simultaneous  Radial  and  Pitch-Angle  Difiiision  and  Dynamic  Outer  Radiation  Belt 
Modeling 

In  the  paper  bearing  this  title  [Chiu  et  al.,  1990b]  we  investigated  (a)  dynamic 
particle  diffusion  processes  caused  by  magnetic  fluctuations  in  the  magnetosphere,  and  (b) 
interacticm  of  radiation-belt  particles  with  the  magnetopause.  We  showed  again  that  the 
general  solution  to  a  simple  form  of  the  bimodal  (radial  and  pitdi-angle)  diffiisirai 
equation,  if  constracted  pieoewise-analytically  with  respect  to  L  value,  can  represent  the 
SCATHA  outer-belt  electron  distribution  with  a  high  degree  of  accuracy  for  up  to  about 
ten  hours  at  a  time.  Representation  of  the  SCATHA  data  in  dus  way  required  us  to  model 
the  diffusion  coefficients  with  adjustable  parameters  whose  numerical  values  are  in  general 
agreement  with  those  obtairwd  in  previous  obsovational  and  theoretical  studies.  Thus,  the 
qjproach  seems  to  form  the  basis  for  the  (fynarruc  modeling  of  hi^-energy  electrons  in  the 
outer  radiation  belt,  at  least  out  to  L  ~  7.  To  elucidate  the  “antnnalous”  behavior  of 
energetic  particles  on  drift  shells  beyond  L  ~  7,  we  performed  single-particle  trajectory 
simulations  (for  protons,  however,  so  as  to  follow  the  particle  gyration  without  usirig 
excessive  amounts  of  computer  time).  From  these  simulations  we  found  several  new 
effects  which  suggest  duu  the  anomalous  behavior  of  electrons  at  L  ^  1  may  have 
involved  encoimters  with  the  magnetopause.  Our  single-particle  trajeaory  simulatirxis  of 
energetic-proton  encounters  with  either  smooth  or  irregular  model  magnetopause  and 
magnetosheath  fields  demcmstrate  the  possibility  of  a  dynamic  recycling  of  radiation-belt 
particles  through  changes  in  the  magnetopause  location  and  through  changes  in  die 
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anq)litudes  of  field  irregulaiities  thereon.  Of  great  interest  in  the  context  of  radiation-belt 
dynamic  modeling  is  that  this  recycling  is  non-difiusive  in  character. 

5.0  Third  Adiabatic  Invariant:  Estimatimi  and  Use 

Organization  of  the  CRRES  outer-belt  electron  data  base  for  modulation  and 
transpon  studies  entails  the  construction  of  time-dependem  profiles  diat  specify  the  {diase- 
space  density  /  as  a  fiinctian  of  the  third  adiabatic  invariant  for  specified  values  of  the 
first  two  adiabatic  invariants  M  and  /. 

A  foimal  confutation  of  the  third  adiabatic  invariant  as  the  amount  of  magnetic 
fiux  enclosed  by  a  particle’s  drift  shell  is  a  time-consuming  procedure  that  seems  too 
expensive  for  routine  use  in  organizing  a  data  base.  Time-saving  alternative  procedures 
developed  in  the  course  of  the  present  work  circumvem  these  difficulties.  One  such 
procedure  entails  line-integration  of  the  vector  potential  on  the  Earth’s  surface  instead  of 
surface-integration  of  the  magnetic  field’s  normal  confonent  over  the  whole  equatorial 
plane.  The  drift  shell  still  needs  to  be  confuted,  but  nufping  it  alcmg  field  lines  to  the 
Eartii’s  surface  (where  the  vector  potential  is  d^le-dominaied,  or  at  least  determined 
essentially  by  currents  inside  the  Earth)  greatly  reduces  the  oonf  lexity  of  the  required 
imegrand. 

The  other  (even  faster)  procedure  sfplied  in  the  present  work  involves  tracing  the 
drift  shell  to  the  nearest  ‘Vefetenoe  longitude”  (typically  either  the  dawn  or  the  dusk 
meridian)  and  determining  the  invariant  label  of  the  correfonding  field  line.  This 
prescr^rm  is  based  on  an  idea  put  forward  by  Roederer  [1970,  p.  107]  that  the  L  value 
of  a  drift  shell  actually  be  defiiKd  as  being  inversely  proportional  to  the  amount  of 
magnetic  flux  enclosed  by  the  shdl.  Magnetospheric  drift  shells  of  very  energetic  particles 
are,  of  course,  essentially  symmetric  across  die  rxxm-midnight  meridional  plane  but  show 
a  inonounced  asymmetry  between  day  and  night  (within  tire  noon-midnight  meridional 
plane  and  elsevdiere).  Moreover,  the  sense  and  degree  of  drift-shell  asymmetry  depends 
on  tire  equatorial  pitch  angle  (or  preferably,  to  put  it  in  adiabatically  invariarrt  terms,  on  the 
value  of  s  vdrere  ihq  is  the  rest  mass).  The  adiabatically  invariartt  K 

parameter  was  introduced  by  Mclhvain  [1966a]  following  a  suggestion  of  Kat^mam 
[1965]. 

The  idea  of  irrvddng  a  ’^reference  longitude”  to  organize  particles  with  respect  to 
drift  shells  seems  to  have  been  introduced  by  Roederer  and  Schuiz  [1969].  However,  the 
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reference-longitude  concept  was  significandy  daiified  by  Schub  [1972,  p.  629]  through  a 
proof  that  the  magnetic  flux  enclosed  by  a  magneticaUy  asymnetic  drift  shell  is  inversely 
proportional  (in  a  first-order  approximation)  to  the  L  value  of  the  field  line  that  intersects 
the  Earth  (or,  to  be  more  precise,  emerges  fiom  the  point  dipole)  at  dawn  or  dusk  along 
the  drift  shell  of  interest.  (The  L  value  of  a  field  line  is  conceived  as  the  limiting  value  of 
r  csc^0  as  r  — »  0,  where  6  denotes  magnetic  colatimde,  so  that  L  is  inversely  proportional 
to  the  Euler  potential  commonly  called  a  or  sometimes  v.  The  Euler  potential  commonly 
called  ^  corresponds  to  the  limiting  value  of  magnetic  longitnde  ^  as  r  — »  0  along  the  field 
line  of  interest.)  Thus,  the  tracing  of  a  drift  shell  to  its  appropriate  reference  longimde 
readily  generates  a  good  approximation  for  the  drird  adiabatic  invariant  <!>. 

To  show  how  the  reference-longitude  ipproach  works  in  practice,  Cladis  et  al. 
[1994]  have  used  it  to  organize  CRRES  energetic-electron  data  fiom  Orbit  270  (executed 
13  November  1990,  ^cfa  was  Day  317  of  the  year).  A  representative  page  of  output 
fiom  the  computer  code  (called  DSTRFN,  for  “distribution  function”)  written  for  this 
purpose  is  shown  in  Table  1.  The  pitch-angle  distribudcHi  ii  diis  case  shows  a  weak 
'‘butterfly”  signature,  widi  an  intensity  maxirmim  at  Oq  »  40°,  even  after  mapping  to  the 
dawn/dusk  reference  longitude.  This  is  a  somewhat  surprising  result,  as  it  seems  to 
indicate  that  the  “butterfly”  signature  is  not  (as  is  widely  believed)  entirely  a  consequence 
of  the  drift-shell  asymmetry.  Data  presorted  by  Pfitzer  et  al.  [1969]  aird  Luhmann  and 
Schulz  [1979]  suggest  a  similar  interpretatior  when  viewed  in  retrospect,  as  the  “butterfly” 
signature  at  0.5-1 .0  MeV  in  those  studies  (based  on  the  same  ATS-1  data  points)  clearly 
extended  about  30°  in  local  time  beyond  the  dawn/dusk  meridian  (i.e.,  to  about  08  hr  and 
16hrLT). 
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Table  1.  Sample  of  output  from  DSTRPN 


YEAR  DAY  OT(SEC)  ORBIT  XKI*  6BO_R(RE)  LAT(DBCI  I.OMC(DEC) 
90  317  68227  270  1  3.488  15.080  80.000 


E(MEV) 

0.760 

0.760 

0.760 

0.760 

0.760 

0.760 

0.760 

0.760 

0.760 

0.760 

0.760 

0.760 

0.760 

0.760 

0.760 


EQ.P.A. (DEG) 
2.0000 
6.0000 
9.6914 
13.5639 
17.4320 
21.2943 
25.1491 
28.9942 
32.8270 
36.6445 
40.4424 
44.2154 
47.9564 
51.6555 
55.2989 


ROEOERERS 

0.0000 

0.0000 

3.7376 

3.7376 

3.7385 

3.7385 

3.7385 

3.7385 

3.7385 

3.7385 

3.7385 

3.7394 

3.7394 

3.7394 

3.7394 


0.760 

58.8666 

3.7394 

0.760 

62.3290 

3.7394 

0.760 

65.6407 

3.7394 

0.760 

68.7309 

3.7394 

0.760 

71.4880 

3.7394 

0.760 

73.7413 

3.7394 

0.760 

75.2577 

3.7394 

0.760 

75.7991 

3.7394 

E(MEV) 

EQ.P.A. (OEG) 

ROEOERERS 

1.156 

2.0000 

0.0000 

1.156 

6.0000 

0.0000 

1.156 

9.6914 

3.7394 

1.156 

13.5639 

3.7394 

1.156 

17.4320 

3.7394 

1.156 

21.2943 

3.7394 

1.156 

25.1491 

3.7394 

1.156 

28.9942 

3.7394 

1.156 

32.8270 

3.7394 

1.156 

36.6445 

3.7394 

1.156 

40.4424 

3.7394 

1.156 

44.2154 

3.7394 

1.156 

47.9564 

3.7394 

1.156 

51.6555 

3.7394 

1.156 

55.2989 

3.7394 

1.156 

58.8666 

3.7394 

1.1S6 

62.3290 

3.7394 

1.156 

65.6407 

3.7394 

1.156 

68.7309 

3.7394 

1.156 

71.4880 

3.7394 

1.156 

73.7413 

3.7394 

1.156 

75.2577 

3.7394 

1.156 

75.7991 

3.7394 

E(MEV) 

EQ.P.A. (OEG) 

ROEOERERS 

1.554 

2.0000 

0.0000 

1.554 

6.0000 

0.0000 

1.554 

9.6914 

3.7394 

1.554 

13.5639 

3.7394 

1.554 

17.4320 

3.7394 

1.554 

21.2943 

3.7394 

1.554 

25.1491 

3.7394 

1.554 

28.9942 

3.7394 

1.554 

32.8270 

3.7394 

1.554 

36.6445 

3.7394 

1.554 

40.4424 

3.7394 

1.554 

44.2154 

3.7394 

K*2(6NEBERS) 

O.OOOE-fOO 

O.OOOE-H)0 

3.802E-02 

1.498E-02 

7.041E-03 

3.663E-03 

2.029E-03 

1.175E-03 

6.984E-04 

4.231E-04 

2.385E-04 

1.591E-04 

9.725E-05 

5.908E-03 

3.333E-05 

2.073E-05 

1.197E-05 

6.640E-06 

3.641E-06 

1.999E-06 

1.134B-06 

7.389E-07 

6.408E-07 

R‘‘2  (CNEBERS) 
0.00OE4-O0 
0.00OE4O0 
3.809BH72 
1.901E-02 
7.046E-03 
3.668E-03 
2.032E-03 
1.176E-03 
6.996E-04 
4.24XE-04 
2.592E-04 
1.591E-04 
9.725E-05 
5.908E-05 
3,333E-05 
2.073E-05 
1.187E-05 
6.640E-06 
3.641E-06 
1.999E-06 
1.134E-06 
7.389E-07 
6.408E-K>7 

K'2 (CNEBERS) 
0.000E400 
O.OOOE-fOO 
3.809E-02 
1.501E-02 
7.046E-03 
3.668E-03 
2.032E-O3 
1.176E-03 
6.996E-04 
4.241E-04 
2.592E-04 
1.591E-04 
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H(Mrr/G) 

F(1/(MEV«NS)«« 

O.OO0E4-O0 

O.OOOE+OO 

0.000E-K>0 

0 . OOOE+00 

5.427E-t-00 

5.039E>01 

1.053£'t-01 

5.270E-01 

1.719E.'01 

5.363E-01 

2.525E401 

5.618E-01 

3.458E-I-01 

5.436E-01 

4.499E-f01 

5.5968-01 

S.627E4-01 

5.5B4E-01 

6.021E-t-01 

5.44BE-01 

8.0S8E-H>1 

5.696E-01 

9.312E-K>1 

5.627E-01 

1.056E+02 

5.578E-01 

1.178E402 

5.440E-01 

1.294E'f02 

5.505E-01 

1.403E-t-02 

5.355E-01 

1.502E-H)2 

5.217E-01 

1. 5898+02 

5.04BE-01 

1.6638+02 

5.022E-01 

1.722E+02 

5.048E-01 

1.765E+02 

5.115E-01 

1.791E+02 

4.953E-01 

1. 8008+02 

5.088E-01 

M(IKV/C) 

F(1/(MEV*NS)»*3 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+00 

O.OOOE+OO 

1.009E+01 

1.3B8E-01 

1.958E+01 

1.443E-01 

3.195E+01 

1.4S1E-01 

4.695E+01 

1.513E-01 

6.429E+01 

1.545E-01 

8.364E+01 

1.530E-01 

1.046E+02 

1.55BE-01 

1.26BE+02 

1.503E-01 

1.498E+02 

1.599E-01 

1.731E+02 

1.588E-01 

1.963E+02 

1.580E-01 

2.190E+02 

1.510E-01 

2.406E+02 

1.527E-01 

2.608E+02 

1.534E-01 

2.792E+02 

1.442E-01 

2.954E+02 

1.444E-01 

3.091E+02 

1.431E-01 

3.201E+02 

1.406E-01 

3.281E+02 

1.459E-01 

3.329E+02 

1.3B2E-01 

3.346E+02 

1.435E-01 

M(ICV/C)  F(1/(MEV*MS)**3 


O.OO0E4-OO 
O.OO0E4-O0 
1.604Ef01 
3.113E401 
S.OBOE'fOl 
7.465E-M)1 
1.022E402 
1.330S4-02 
1 . 6C3E-I-02 
2. 01 6E  44)2 
2.382E4-02 
2 . 7S3E402 


O.OOOE-t-00 
O.OOOE'fOO 
4.202E-02 
4.431E-02 
4  .265E-02 
4 .490E-02 
4.557E-02 
4 .64SE-02 
4.869E-02 
4.020E-02 
5.064E-02 
S.013E-02 
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ABSTRACT 


Intensities  of  geomagnetically  trapped  outer-zone  radiation-belt  electrons  typically 
decrease  (often  by  orders  of  magnitude)  as  the  main  phase  of  a  magnetic  storm  develops, 
then  increase  at  least  back  to  pre-storm  levels  (and  often  higher)  during  recovery  phase. 
The  “cleanest”  events  are  those  in  which  a  single  decrease  in  (from  values  near  zero 
toward  values  around  -100  nT)  is  followed  by  a  smooth  recovery  without  further 
fluctuation.  The  present  work,  a  systematic  study  of  the  ^ical  relativistic-electron 
response  to  such  changes  in  Z)^^,  invokes  a  CRRES  data  base  including  about  IS  of  such 
“clean”  events  (as  well  as  some  less  “clean”  events)  that  occurred  during  the  lifetime  of 
CRRES.  The  purpose  of  this  study  is  to  (fetermine  the  extent  to  which  the  ^ical  outer- 
belt  electron  response  to  D^i  is  essentially  adiabatic  (consistent  with  conservation  of  the 
three  adiabatic  invariants  during  a  global  inflation  and  subsequent  relaxation  of  drift  shells) 
and  conversely,  the  residual  extent  to  which  the  typical  starmtime  decrease  in  energetic- 
electron  intensities  involves  a  genuine  loss  of  trapped  particles.  Two  approaches  to  the 
analysis  of  relevant  CRRES  data  arc  attempted  in  the  present  study.  The  more  traditional 
approach  makes  use  of  a  model  for  the  magnetic-field  perturbation  AB  produced  by  the 
ring  current.  The  novel  feature  of  our  tiKxlel  is  that  both  the  equatorial  AB(r)  and  the 
corresponding  magnetic  flux  enclosed  by  a  drift  shell  of  radius  r  are  expressible 
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analytically  in  teims  of  the  model  parameters.  The  less  traditional  approach  makes  use  of 

the  CRRES  energetic-electron  and  magnetic-field  data  so  as  to  reveal  what  radial 
variation,  for  the  ratio  (d  In  B^/d  In  between  the  fiactional  changes  in  B  seen  by  the 

particle  and  by  an  in  situ  observer,  respectively,  would  be  required  in  order  that  the 
observed  electron-flux  modulation  might  be  purely  adiabatic  in  origin.  The  purpose  of 
having  tried  this  alternative  scheme  was  to  learn  whether  it  might  offer  further  insight  into 
stormtinie  responses  of  energetic  electrons  to  variations,  and  whether  such  stormtime 
responses  might  help  to  reveal  the  spatial  structure  of  the  magnetic-field  perturbation  AB 
produced  by  the  actual  ring  current  The  inferred  (d  In  B^/d  In  B^)  often  closely  resembled 

that  derived  from  the  radial  variation  of  AB  in  standard  ring-current  models,  as 
surmnarized  by  our  generic  analytical  representation,  but  in  most  such  cases  the  AB(r) 
actually  observed  was  not  in  accord  with  the  standard  truxiels.  The  more  traditional 
analysis  of  the  stormtime  energetic-eiectron  response  to  changes  in  Z),,,  based  on  the 
adiabatic  transformation  of  radial  profiles  and  energy  spectra  observed  between  storms, 
suggests  that  adiabatic  response  to  a  realistically  modeled  AB(r)  can  account  for  anywhere 
from  zero  to  100%  (but  typically  less  than  half)  of  the  observed  stormtime  decrease  in 
logarithmic  electron  flux  as  decreases  toward  mote  negative  values,  but  scatter  among 
the  data  points  mates  it  difficult  to  mate  a  more  quantitative  estimate  than  this. 
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INTRODUCnON 

A  major  challenge  in  the  observational  study  of  panicle  data  from  the  outer 
radiation  belt,  and  in  the  construction  of  reliable  time-dependent  models  of  the  radiation 
environment,  is  the  separation  of  adiabatic  (reversible)  variations  horn  transpon-rtiaed 
(diffusive)  phenomena.  This  challenge  is  especially  foimidable  during  geomagnetic 
storms,  which  are  characterized  on  the  one  hand  by  unusually  strong  radial  diffusion  (a 
transport  process)  and  on  the  other  hand  by  an  unusually  strong  ring  current  (which  is 
responsible  for  major  adiabatic  variations  in  energetic-paiticle  intensities).  Indeed,  radial 
diffusion  of  plasma-sheet  and  ring-current  ions  is  an  important  process  for  producing  the 
stormtime  ring-current,  to  wdiich  the  higher-energy  radiation-belt  particles  respond  quasi- 
adiabatically.  Although  the  radiation-belt  electrons  of  interest  in  the  present  smdy  are 
much  more  energetic  than  typical  ring-current  particles,  they  are  (of  course)  subject  to 
essentially  the  same  forces  and  drifts. 

We  focus  in  the  present  work  on  the  adiabatic  responses  of  radiation-belt  electron 
intensities  to  variations  in  the  geomagnetic  ring-current  index  We  characterize 
adiabatic  responses  here  as  “shon-term”  aixl  diffusive  responses  as  “long-term.”  This 
distinction  is  based  on  the  premise  that  adiabatic  response  is  conceived  as  essentially 
immediate,  whereas  diffusive  response  is  conceived  as  a  cumulative  fdicnomenon. 
However,  it  is  clear  that  the  same  magnetospheric  event  (e.g.,  a  sudden  commencement) 
could  easily  entail  elements  of  both  processes.  The  response  of  radiation-belt  electrons  to 
the  changes  in  D„  that  draracteristically  follow  a  sudden  commencement  is  more  likely  to 
be  adiabatic,  since  (a)  the  time  scale  on  £>^  varies  is  mote  likely  to  exceed  the  drift 
periods  of  the  particles  of  interest  and  (b)  the  magnetic  field  produced  by  the  ring  current 
is  more  nearly  symmetric  in  azimuth  than  is  the  magnetic  disturbance  associated  with  a 
sudden  commencement. 
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In  any  case,  the  standard  Idneniatical  framewoik  for  radiation-belt  modelling  is 
provided  by  the  three  adiabatic  invariants  (Af,  J,  <i>),  vdiich  arc  proportional  (respectively) 
to  canonical  action  integrals.  These,  together  with  their  conjugate  i^iases  (which  are  often 
suppressed  or  averaged-over),  constitute  the  phase  space  occupied  by  the  electron 
population.  A  con^lete  description  of  the  electron  radiation  belts  is  achieved  in  piinc^le 
by  specifying  the  phase-space  density  as  a  function  of  time.  When  the  magnetosphere  is 
disturbed,  as  during  the  the  onset  and  main  phase  of  a  geomagnetic  storm,  the  three 
canonical  action  variables  are  unlikBly  to  remain  constant.  (They  are,  after  all,  only 
“adiabatically”  invariant.)  The  adiabatic  invariants  remain  largely  valid,  however,  as 
kinematical  variables  with  respect  to  whidi  the  evolution  of  radiation-belt  electron  j^iase- 
space  density  is  most  easily  described.  Only  by  organizing  particle-flux  data  with  respect 
to  the  three  adiabatic  invariants  can  we  determine  the  extent  to  which  non-adiabatic 
transport  has  affected  the  panicle  distiibuticm.  Thus,  the  adiabatic  invariants  constimte  an 
essential  framework  for  any  dynamical  study  of  the  radiation  belts. 

The  level  of  geomagnetic  activity  aiKl  the  rate  of  energetic-patticle  transpon  are 
both  characterized  [Mayaud,  1980]  by  the  quasi-logarithmic  index  Kp,  >^ch  is  computed 
separately  for  each  3-hr  time  interval  and  thus  eight  times  per  day.  The  sum  of  the  eight 
Kp  values  for  a  given  day  (UT)  is  denoted  ZKp.  Those  days  for  ^^ch  IXp  is  lowest  and 
highest  are  designated,  respectively,  the  most  quiet  ((JQ)  and  most  disturbed  (DD)  days  of 
the  month.  The  strength  of  a  geomagnetic  storm  and  the  distortion  of  the  geomagnetic 
field  during  it  are  better  indicated,  however,  by  the  ring-current  index  This  is 
compiled  at  hourly  intervals  from  data  provided  by  a  network  of  several  low-latitude 
magnetic  observatories.  The  Dg^  index  is  a  direct  measure  of  the  horizontal  (H)  magnetic- 
field  perturbation  (relative  to  a  nominally  quiescent  baseline)  imposed  at  the  average  near- 
equatorial  ground  station  by  the  ring  current  and  by  telluric  currents  tiiereby  induced. 
Values  of  order  -200  nT  to  -250  nT  for  are  characteristic  of  major  geomagnetic 
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storms.  An  even  more  negative  value  (D,t  =  -298  nT)  was  seen  during  the  huge  storm  of 
24-25  March  1991. 

Koons  and  Gomey  [1991]  have  reported  what  amounts  to  an  inverse  conelation 
between  relativistic-electron  fluxes  at  synchronous  altitude  and  the  weighted  sum  of  the 
Kp  index  over  the  previous  several  days.  A  plausible  interpretation  of  this  finding  is  that 
Kp  should  be  a  leading  indicator  (although  not  by  as  much  as  several  days)  for  D^^. 
(Enhanced  Kp  accompanies  enhanced  transport  of  20-200  keV  ions  into  the  ring-cunent 
region,  where  they  generate  a  negative  However,  an  increasingly  negative  D,^  should 
decrease  the  energies  of  already-trapped  relativistic  elections  and  thereby  reduce  the 
count  rates  of  such  electrons  in  orbiting  particle  detectors  (which  have  fixed  energy 
thresholds,  of  course).  Decay  of  the  stormtime  ring  current  during  recovery  phase  should 
restore  relativistic  electron  count-rates  to  at  least  their  pre-storm  values  (and  perhaps 
above,  in  view  of  the  stonntime  transport  that  such  electrons  would  have  experienced  in 
response  to  the  sudden  commencement). 

In  the  present  woric  we  measure  the  extent  to  which  adiabatic  response  to  temporal 
variations  in  can  account  quantitatively  for  temporal  variations  in  relativistic-electron 
intensities  observed  in  the  radiation  belts  during  the  SPACERAD  (space-radiation)  phase 
of  the  USAF/NASA  CRRES  mission  (Combined  Release  and  Radiation  Effects  Satellite) 
firom  September  1990  through  October  1991.  The  prototype  for  our  study  is  found  in  the 
work  of  Soraas  and  Davis  [1968],  who  organized  radiation-belt  proton  observations 
[Davis  and  Williamson,  1963]  into  phase-space  density  profiles  by  using  a  sinqtle  model 
for  the  magnetic-field  distortion  associated  with  Related  studies  pertaining  to  ring- 
current  ions  have  been  performed  by  Lyons  and  Williams  [1976]  and  followed-up  by 
Williams  [1981],  but  the  basic  idea  (adiabatic  modulation  of  energetic-particle  fluxes  by 
Djt)  seems  to  have  originated  with  Dessler  and  Karplus  [1961]  and  was  pursued  also  by 
Mcllwain  [1966h]. 
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The  fact  that  relativistic  electron  fluxes  in  the  outer  radiation  belt  undergo  drastic 
variations  (often  amounting  to  several  orders  of  magnitude)  in  response  to  geomagnetic 
activity  has  been  known  since  the  early  days  of  outer-belt  observation  [e.g..  Brown  et  al., 
1968].  Although  outer-belt  electron  intensities  vary  much  more  strongly  in  response  to 
geomagnetic  activity  than  do  (for  example)  inner-zone  trapped-proton  intensities,  a 
physical  description  of  the  dynamics  involved  must  be  developed  along  similar  lines:  use 
of  the  three  adiabatic  invariants  (//)  as  kinematic  variables,  specification  of  charged- 
particle  distributions  (phase-space  densities)  as  time-dependent  functions  of  the  {/,},  and 
characterization  of  particle  transport  in  terms  of  the  violation  of  one  or  more  of  the 
adiabatic  invariants.  (Diffusion  and  “dynamical  fiiction”  in  phase  space  encompass 
transport  in  ordinary  space  as  well  as  in  energy  and  pitch  angle.)  As  any  such  transport  is 
superimposed  on  a  kinematical  firameworic  that  treats  the  adiabatic  invariants  as  if  they 
were  conserved,  this  framework  serves  to  order  and  organize  the  phase-space  density 
profiles  of  outer-belt  electrons  even  if  some  of  the  adiabatic  invariants  are  weakly  violated 
(such  as  during  the  course  of  a  geomagnetic  storm).  Thus,  it  makes  sense  to  model  the 
phase-space  density  as  a  function  of  the  three  adiabatic  invariants,  even  through  an  event 
during  which  the  invariants  are  not  strictly  conserved.  {Lanzerotti  et  al.  [1970]  organized 
some  outer-belt  electron  data  into  phase-space  density  profiles  for  the  purpose  of 
extracting  transport  coefficients  but  treated  the  magnetic  field  as  static  in  time  and  thus  did 
not  invoke  as  a  modeling  parameter.) 

Given  sufficient  coverage  of  phase  space  with  an  ^ropriate  but  time-limited  data 
set,  one  can  successfully  focus  upon  anecdotal  aspects  of  outer-belt  dynamics,  such  as 
radial  diffusion  and/or  pitch-angle  diffusion  on  a  time  scale  of  S-20  days  [e.g.,  Lanzerotti 
et  al.,  1970;  West  et  al.,  1981;  Chiu  et  al.,  1988,  1990].  Such  episodic  interpretations  of 
partial  data  sets  are  quite  valuable.  This  is  the  perspective  adopted  in  the  present  work. 


29 


7 


although  our  long-range  purpose  is  to  model  the  variation  of  radiation-belt  particle 
intensity  as  a  continuous  function  of  the  spatial  coordinates,  energy,  pitch-angle,  and  time. 

In  order  to  achieve  the  more  complete  dynamical  nxxlel  of  the  outer-belt  electron 
environment  and  its  response  to  geomagnetic  activity,  one  must  (of  course)  look  for  a 
long-term  description  that  takes  account  of  conservative  as  well  as  non-conservative 
processes.  Since  the  geomagnetic  field  undergoes  drastic  variations  with  geomagnetic 
activity  in  the  outer-belt  region,  a  complete  description  of  the  outer-belt  response  must 
involve  both  electron-flux  measurements  and  corresponding  geomagnetic  field 
measurements  on  all  the  time  scales  contained  in  the  data  set  Soraas  and  Davis  [1968] 
have  shown,  in  a  study  involving  outer-zone  protons,  how  to  achieve  such  a  description  in 
principle.  We  show  here  that  a  similarly  systematic  ordering  of  the  outer-zone  electron 
environment  could  be  achieved  by  means  of  our  CRRES/SEP  data  base. 

DATA  BASE 

The  CRRES  satellite  was  launclvKl  25  July  (Day  206)  1990  into  a  highly  elliptical 
orbit  (perigee  =  350  km,  apogee  =  33584  km  =  5.27  that  reached  almost  to 
geosynchronous  altitude  (=  5.59  and  provided  useful  data  until  12  October  (Day  285) 
1991  [Johnson  and  Kierein,  1992;  Gussenhoven  and  Mullen,  1993].  [Collections  of 
CRRES  papers  have  been  published  in  IEEE  Trans.  Nucl.  Sci.,  38  (6),  December  1991;  in 
/.  Spacecrerft  Rockets,  29  (4),  July-August  1992;  and  in  IEEE  Trans.  Nucl.  Sci.,  40  (2), 
April  1993].  The  plasma-release  component  of  the  CRRES  mission  involved  an  array  of 
24  chemical  canisters  expended  during  the  course  of  three  separate  campaigns.  These 
releases  are  not  of  interest  in  the  present  context.  The  present  work  is  related  instead  to 
the  SPACERAD  component  of  the  CRRES  mission,  which  had  the  purpose  of  making 
detailed  observations  on  the  Earth’s  radiation  belts.  Indeed,  one  of  the  major  scientific 
objectives  of  CRRES  was  to  elucidate  the  dynamics  of  the  high-energy  electrons  in  the 
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outer  radiation  belt  by  means  of  a  corrplement  of  modem  instruments  [Vampola,  1992]. 
Here  we  report  on  temporal  variations  of  relativistic-clectron  intensities  (£  ^  300  keV) 
observed  at  L  ~  4-7  in  the  outer  radiation  belt  by  the  ONR-307  Spectrometer  for 
Electrons  and  Protons  (SEP)  aboard  the  CRRES  satellite  [Nightingale  et  ai,  1992].  (The 
“L”  value  invoked  in  this  paper  is  a  nominal  one,  defined  below  for  the  purpose  of 
reference,  but  not  used  as  a  physical  parameter  in  actual  calculations.) 

The  data  base  for  the  present  study  has  three  components:  (1)  relativistic  electron 
fluxes  fi^om  the  aforementioned  CRRES  ONR-307-3  SEP  instrument,  (2)  the  in  situ 
magnetic  field  vector  provided  by  the  onboard  CRRES  magnetometer,  and  (3)  the 
geomagnetic  index  D^^. 

Our  present  study  is  event-oriented.  It  comprises  a  systematic  analysis  of  the 
typical  relativistic-electron  response  to  changes  in  associated  with  magnetic  storms  of 
various  intensities.  The  “cleanest”  events  arc  those  in  which  a  single  decrease  in  (from 
values  near  zero  toward  values  around  -100  nT)  is  followed  by  a  smooth  recovery 
without  further  fluctuation.  The  data  base  for  the  present  work  includes  (see  Table  1)  all 
13  of  such  “clean”  events  that  occurred  during  the  14-month  lifetime  of  CRRES,  plus  a 
few  major  storms  that  were  less  “clean”  in  this  respect.  Separate  columns  in  Table  1 
indicate  the  most  negative  value  of  Dgf  attained  during  each  storm  in  the  data  base,  as  well 
as  the  factor  by  which  electron  fluxes  in  the  561-958  keV  energy  channel  (nominal  energy 
670  keV,  as  is  explained  below)  were  typically  reduced  during  the  conesponding  storm. 

The  CRRES  ONR-307-3  SEP  instrument  [Nightingale  et  al.,  1992]  is  similar  to 
the  energetic-partcle  detector  flown  on  the  1979  SCATHA  spacecraft  [Reagan  et  al., 
1981].  It  consists  of  a  stack  of  three  identical  sensors  (called  A,  B,  and  C).  Each  has  a  3° 
angular  field  of  view  (FWHM).  The  three  sensors  were  oriented  at  40°,  60°,  and  80° 
(respectively)  relative  to  the  spacecraft  spin  axis,  and  each  covered  the  energy  range  42 
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keV  to  5  MeV  in  12  energy  channels.  Table  2,  taken  from  Nightingale  et  al.  [1992], 
shows  the  modes  of  operation  and  enei^  coverage  for  the  SEP  instrument  The  electron 
data  base  for  the  present  study  comes  from  the  pitch-angle  and  energy  distributions 
provided  by  the  various  electron  modes  of  the  SEP  instrument  Energy  spectra  deduced 
from  the  present  data  base  have  been  compared  extensively  with  data  from  the  Nfedium 
Energy  Analyzer  (MEA)  instrument  on  CRRES  [Vampola  etal.,  1992]  and  found  to  be  in 
good  agreement  with  the  latter  except  during  major  storms,  such  as  the  hist(»ic  storm  of 
24  March  1991,  when  our  SEP  instrument  seemed  to  have  been  affected  significantly  by 
bremsstrahlung  contamination.  Accordingly,  we  restrict  oar  primary  attention  here  to 
events  for  which  |f>^|  ^150nT. 

Together  with  similar  CRRES  instruments  [e.g.,  Vampola  et  al.^  1992],  our  data 
set  based  on  the  ONR-307/SEP  instrunoent  [Nightingale  et  al.,  1992]  coir^lements 
previous  near-geosynchronous  satellite  data  sets  of  much  longer  duration  [e.g..  Baker  et 
al.,  1978,  1990;  Reagan  et  al.,  1987;  Chiu  et  al.,  1988;  1990]  and  offers  a  qualitative 
improvement  upon  particle  data  with  more  limhed  energy  nd/or  pitch-angle  resolution 
[e.g..  Brown  et  al.,  1968;  West  et  al.,  1981;  Sibeck  et  al.,  1987]. 

We  have  obtained  CRRES  magnetometer  data  [Singer  et  al.,  1992]  with  verified 
calibration  for  this  study  from  agency  tapes  supplied  through  the  CRRES  data-distribution 
network.  The  stormtime  scalar  magnetic-field  intensities  B  (s  |B|)  used  in  the  present 
study  are  from  actual  magnetometer  measurements  rather  than  from  models.  The 
stormtime  ring  current  (as  measured  by  the  geomagnetic  index  can  have  a  strong 
influence  on  B  in  the  region  of  the  outer  radiation  belt,  but  the  most  widely  available  time- 
dependent  magnetosphoic  B-field  model  [Tsygmenko,  1989]  expresses  time-dependence 
through  Kp  instead  of  through  D^^.  The  B-field  model  provided  by  the  CRRES  data 
network  is  the  most  recent  Olson-Pfitzer  model,  which  is  an  improved  version  of  the 
model  described  by  Olson  and  Pfitzer  [1974].  This  is  always  in  good  agreement  with  the 
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GIRES  magnetometer  data  at  low  altitudes,  where  the  Earth’s  main  (internal)  field 
dominates.  However,  the  Olson-Pfitzer  model  pertains  to  time-averaged  geomagnetic 
conditions  and  so  contains  no  time-dependence  corresponding  to  that  of  D^^. 

We  use  Djt  and  not  Kp  as  our  reference  index  for  identifying  events  of  interest  to 
the  present  study.  As  this  is  a  study  of  trapped-partick  responses  to  geomagnetic 
variability,  we  need  a  geomagnetic  index  that  specifies  the  present  geomagnetic  field 
configuration.  Having  made  correlations  of  CRRES/SEP  relativistic  outer-belt  electron 
fluxes  (£  ^  300  keV)  with  the  Kp  and  D^t  iniices  over  the  lifetime  of  (ZRRES,  we  found 
(as  others  had  found  previously)  only  a  vague  correlation  of  these  fluxes  with  Kp  but  a 
very  clear  connection  between  relativistic  electron  intensities  and  Z)^^.  (This  is  understood 
to  mean  that  Kp  indicates  only  the  present  rate  of  transport  for  both  ring-current  and 
radiation-belt  particles,  whereas  £),(  indicates  the  present  energy  content  of  the  ring 
current  and  thus  the  present  configuration  of  the  magnetospheric  B  field.) 

Figure  1  shows  a  representative  50-day  coirqianison  of  energetic-electron 
intensities  (.J^)  in  two  energy  channels  (E  =  561-958  keV  and  £  -  958-1355  keV)  for  4.5 

^  L  ^  6.5,  together  with  the  geomagnetic  index  (available  as  a  histogram  with  1-hr 
time-resolution)  for  the  same  50-day  interval  For  typical  energy  spectra  (see  below) 
encountered  during  the  CRRES  mission,  an  energy  of  670  keV  was  representative  of  the 
561-958  keV  channel  and  an  energy  of  1090  keV  was  representative  of  the  958-1355 
keV  channel 


DATA  ORGANIZATION  AND  INTERPRETATION 

The  CRRES  relativistic-electron  data  base  spans  a  considerable  range  of  magnetic 
latitudes.  This  means  that  the  SEP  pitch-angle  distributions  were  obtained  at 
measurement  points  corresponding  to  various  locations  along  the  field  lines  of  interest. 
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To  order  the  data  base  so  that  all  the  locally  measured  pitch  angle  distributions  have  a 
uniform  reference,  we  “map”  all  the  in  situ  pitch-angle  distributions  to  the  magnetic 
equator  (locus  of  minima  in  B  =  I B  I  along  field  lines).  Some  typical  results  are  shown  in 
Figure  2,  where  plus  signs  (+)  correspond  to  local  pitch  angles  and  boxes  (□)  to 
equatorial  pitch  angles.  Figure  3  shows  an  example  of  a  “butterfly”  pitch-angle 
distribution  [cf.  West  et  al.,  1973],  so-called  because  of  its  four-lobed  appearance  on  a 
360°  polar  plot.  The  CRRES  data  network  provides  a  best-fitting  magnetic  field 
according  to  the  Olson-Pfitzcr  model,  and  we  use  this  for  any  required  mappings.  The 
Olson-Pfitzer  model  always  agrees  well  with  CRRES  magnetometer  data  at  the  lower 
altitudes,  where  the  Earth’s  main  (i.e.,  internal)  field  dominates.  The  model  is  quite 
reliable  also  in  the  region  occupied  by  the  outer  radiation  belt  during  geomagnetically 
quiet  time  intervals. 

In  this  paper  we  obtain  “equatorial”  pitch-angle  distributions  finom  off-equatorial 
measurements  by  tracing  each  field  line  of  interest  (a  curve  everywhere  tangent  to  the 
direction  of  B)  back  to  the  magnetic  equator  (locus  of  minima  in  B  along  field  lines)  in  the 
best-fitting  Olson-Pfitzer  field  model.  This  nuqpping  reveals  a  value  for  Bq,  defined  as  the 
minimum  value  of  B  (s  |  B  { )  along  the  field  line  of  interest,  and  thus  allows  us  to  make 
meaningful  comparisons  among  pitch-angle  distributions  obtained  at  different  magnetic 
latitudes  as  CRRES  cuts  through  the  outer-belt  region.  We  use  the  nominal  Olson-Pfitzer 
“L-value”  of  a  pitch  angle  distribution,  however,  only  as  a  reference  parameter  for 
comparing  local  pitch-angle  distributions  measured  at  different  places.  We  do  not  use  this 
“L-value”  for  dynamical  calculations  here.  We  also  do  not  use  the  static  Olson-Pfitzer 
model  field  for  dynamical  calculations;  we  use  only  measured  magnetic-field  values  and 
our  model  ring-current  field  (see  Appendix)  for  such  purposes.  All  our  spatial  references 
to  CRRES  equatorial  crossings  are  given  in  Earth  radii 
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PHASE-SPACE  DENSITY 

The  phase-space  density  /  is  well  known  to  be  equal  to  the  ratio  of  differential 
unidirectional  flux  JJJE)  to  p^,  where  p  is  the  particle  nnomentum  and  a  s  cos'Kp*  B )  is 

the  local  pitch  angle.  Uouville’s  theorem  asserts  that /is  a  conserved  quantity  in  the  sense 

that  d//d/  s  0  along  a  phase-space  trajectory  consistent  with  Hamiltonian  mechanics. 
Thus,  the  mirroring  particle  flux  is  given  by 

yj.  =  2/noAfB/,  (1) 

where  A  is  the  local  magnetic-field  intensity,  ilf  is  the  first  adiabatic  invariant,  and  /hq  is  the 
rest  mass.  Thus,  any  process  (consistent  with  Hamiltonian  mechanics)  which  causes  the 

mirror-point  field  of  a  particle  to  vary  with  time  will  at  first  sight  necessarily  cause  the 
mirror-point  value  of  to  vary  as  B„  at  fixed  Af.  However,  it  is  customary  (because  of 

the  way  charged-particle  instruments  ate  necessarily  designed)  to  regard  as  a  function 

of  particle  energy  (£),  equatorial  pitch  angle  (Oq),  and  L  value.  It  follows  from  (1)  that 
the  observed  at  a  given  energy  E  (and  thus  at  a  given  value  of  the  product  A/B„) 

should  vary  approximately  as  B^  with 

<T  =  -(dlny?ainA#)jfj,,  (2) 

where  K  (s  is  an  energy-independent  quantity  derived  [Kai^nann,  1965; 

Mcllwain,  1966a;  Roederer,  1970]  from  the  first  two  adiabatic  invariants  {M  and  J). 
However,  the  B  with  which  the  observed  JJ.E1  should  thus  vary  is  the  value  of  B  (=  B^)  at 

the  guiding  center  of  the  mirroring  particle.  This  is  not  necessarily  equal  to  the  value  of  B 
(=  B^)  "observed”  at  CRRES,  since  the  indutxd  electric  field  associated  with  a  time- 

varying  B  would  have  moved  the  particle’s  guiding-center. 
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Upon  encounteiing  a  similar  ambiguily  while  trying  to  inteipret  some  ATS-1  data, 
Barfield  et  al.  [1971]  deduced  that  the  relative  logarithmic  modulation  (particle  flux  vs 
magnetic  intensity)  observed  at  the  spacecraft  should  be  given  by 

(A  In  JJA  In  =  { 1  -  [(7  +  l)/27j(3  In  JJd  In  £)icxKd  In  5^/d  In 

+  0  In  Jjd  In  [1  -  (d  In  B^J(i  In  .  (3) 

where  7  is  the  ratio  of  relativistic  mass  (m)  to  rest  mass  (mo).  The  logic  behind  (3)  is  that 
the  observed  7^  ^  ^  ^  pre-storm  phase-space  density  but 

evaluated  at  a  different  M  (and  possibly  at  a  diffisrent  L)  confonning  to  the  desired  energy 
at  the  spacecraft  and  to  the  geometry  of  the  field  perturbation.  In  a  modulation  that 

hypothetically  moved  the  particles  and  conqnessed/rarefied  the  field  azimuthally,  for 
exanq>le,  the  change  in  B  seen  at  the  spacecraft  (ABq)  would  be  equal  to  the  change  in  B 

experienced  by  the  particle  at  its  guiding  center  (Afi^).  In  this  case  the  pre-storm  phase- 

space  density  of  interest  would  be  that  corresponding  to  a  different  Af  (such  that  A  In  Af  = 
-  A  In  fl^)  at  the  same  L  value,  and  the  value  of  that  phase-space  density  would  be  such 

that 

Ain/  =  (dln//AlnAf)jcx(^^^  =  [1  -  (dln/x/9lnP^)js:xl(AlnB^) .  (4) 

The  factor  [(7  +  l)/2y]  multiplying  0  In  JJB  In  £)^x  “  (3)  represents  (d  In  E/d  In  p^), 
whereas  the  factor  O  In  JJd  In  E)i^j^  represents  the  observed  quiet-time  errergy  spectrum 
and  =  2m^B^, 

The  first  term  on  the  right-hand  side  of  (3)  is  thus  reasonably  self-explanatory  in 
terms  of  (1)  and  (2).  The  second  term  on  the  right-hand  side  of  (3)  contributes  whenever 
temporal  variation  of  B  induces  drift  shells  to  move  “radially”  past  the  spacecraft,  in  which 
case  the  observer  necessarily  samples  the  underlying  “radial  gradient”  (9  In  JJd  In  B^i^^ 

of  the  particle  distribution  as  well  as  the  eirergy  spcctram.  In  this  case  the  change  in  B 
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experienced  by  the  particle  typically  differs  from  the  change  in  B  seen  at  the  spacecraft,  so 
that  (d  In  B^/d  In  #  1  in  (3).  For  example,  a  siiiq>lifiBd  model  in  which  B  is  weakly 
perturbed  by  adding  a  uniform  but  time-dependent  AB  =  parallel  (or  antiparallel)  to 
the  dipole  axis  is  well  known  (and  easily  shown:  see  Appendix)  to  yield  (d  In  In  B^)  = 

S/2  for  equatorially  mirroring  particles  (Le.,  for  K  =■  0).  Interpretation  of  particle-flux 
modulation  as  a  largely  adiabatic  consequence  of  variations  in  ring-current  intensity  is  thus 
inherently  model-dependent  on  account  of  tire  factor  (d  In  B^/d  In  B^).  Estimates  of  the 

quantity  (d  In  B^/d  In  B^)  for  a  reasonable  model  of  the  magnetic  fidd  generated  by  the 
ring  current  are  developed  in  the  Appendix  for  future  use,  as  well  as  for  comparison  with 
results  obtained  here. 

Since  the  ATS-1  satellite  which  had  provided  the  data  analyzed  by  Barfield  et  al. 

[1971]  was  in  geosynchronous  orbit,  there  was  littie  ambiguity  as  to  meaning  of  the 
quantity  (A  In  JJti  In  B^j^^  whose  expected  value  is  specified  by  (3):  Instruments 
onboard  ATS-1  had  provided  direct  measurements  of  both  B^  and  as  functions  of  time, 

and  the  observed  modulation  of  each  showed  a  period  -- 100  sec  at  what  was  essentially  a 
fixed  spatial  location.  Both  drere  and  here,  the  observed  nnodnlation  of  is  regarded  as  a 

consequence  of  dynamical  processes  that  had  led  to  the  observed  modulation  of  B^. 

However,  interpretation  of  the  CRRES  data  by  rteans  of  (3)  is  conq>licated  by  the  fact 
that  CRRES  has  a  highly  ellqrtical  orbit  whose  period  is  comparable  to  the  time-scale  of 

the  physical  process  under  consideration  (growth  and  decay  of  the  stormtime  ring 
current).  Thus,  we  need  to  sort  out  the  spatial  variations  (rf  B^  and  (which  CRRES 

would  have  encountered  as  terrqroral  variations  even  in  a  static  magnetosphere)  from  truly 
in  situ  temporal  variations  of  B^  and  (^^riuch  would  have  resulted  from  growth  and 
decay  of  the  stormtime  ring  current).  Accordingly,  we  deduce  A  In  B^  and  A  In  here  by 
comparing  measurements  of  B^  and  of  Jj_  made  on  orbits  at  most  a  few  days  apart, 

sampling  as  nearly  as  possible  the  same  spatial  location  for  eadi  comparison. 
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RADIAL  PROFILES  AND  ENERGY  SPECTRA 

As  quict-timc  values  for  (9  In  JJd  In  and  (9  In  JJd  In  needed  in 

order  to  test  (3)  against  our  stormtime  CRRES  data,  we  have  made  suitable  fits  to  the 
radial  profiles  (at  fixed  energy)  and  energy  spectra  (at  various  L  values)  observed  on  orbits 
executed  during  the  relatively  quiet  time  intervals  before  the  respective  storms  of  interest 
For  our  first  example  we  decided  to  focus  on  the  moderate  storm  of  27-28  November 
1990  (Days  331-332  of  the  year).  A  representative  radial  profile  (observed  vs 

observed  B)  fi'om  Orbit  297  (executed  1946-0537  UT,  24—25  November  1990,  Days 
328-329  of  the  year)  are  shown  in  Figure  4.  (CRRES  orbits  begin  and  end  at  perigee  by 
convention  and  are  numbered  accordingly.  By  “observed”  B,  we  mean  the  Pythagorean 
sum  of  the  three  spacecraft-oriented  vector  corrqxinents  of  B  reported  by  the  onboard 
CRRES  magnetometer.)  Geomagnetic  conditions  were  very  quiet  during  Orbit  297,  and 
the  value  of  Id^J  remained  <  10  nT  for  the  whole  day.  The  outbound  pass  found 
CRRES  at  0.7®  to  1.5®  nominal  (dipolar)  magnetic  latitude  \  for  4.5  <  L  <  5.5;  the 
inbound  pass  found  CRRES  at  X,  -  5.0®  to  5.8®  over  the  same  range  of  L  values.  As  the 
difference  in  magnetic  latitude  between  inbound  and  outbound  encounters  with  the  same  L 
value  may  have  been  significant  for  our  conterrqilated  analysis,  we  made  separate  (as  well 
as  joint)  least-squares  cubic  fits  (logjo/j.  vs  logio®)  for  the  inbound  and  outbound 

segments  of  the  same  orbit  Results  for  the  £  =  561-958  keV  energy  channel  (identified 
by  the  nominal  energy  E  -  670  keV)  are  shown  in  Figure  4.  Since  the  observations  were 
nearly  equatorial,  the  displayed  range  of  B  (100—1000  nT)  corresponds  to  6.7  ^  L  ^  3.1, 
the  maximum  radiation  intensity  having  been  attained  at  B  =  390  nT  (L  «  4.3).  We  use  the 
joint  fit  (solid  curve,  inbound  and  outbound  data  lumped  together)  for  interpreting  data 
fi'om  the  27-28  November  1990  storm  in  the  present  work. 
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Representative  energy  spectra  at  L  »  4.5  and  L  »  5.5  are  shown  in  Figure  5.  We 
have  fitted  such  spectra  widi  exponential  functions  of  the  farm 

J^{E)  =  Cexp(-E/Eo)  (5) 

at  selected  L  values  and  have  thus  determined  corresponding  values  for  Eq.  Some 
representative  results  for  Eq  are  plotted  iri  Figure  6.  The  superimposed  parabola  suggests 
that  Eq  can  be  fitted  successfully  by  a  quadratic  function  of  L.  (The  scatter  in  values  for 
Eq  at  L  ^  6.5  is  unin^rtant  for  applications  treated  here.) 

Our  functional  fits  to  quiet-time  radial  profiles  and  energy  spectra  could  be  used  to 
constma  phase-space  density  profiles  [e.g.,  Lamerotti  et  al.,  1970]  at  fixed  M  and  J  (first 

two  adiabatic  invariants).  Here,  however,  we  just  differentiate  the  functional  fits 
analytically  so  as  generate  numerical  values  of  (d  In  Jjd  In  and  (d  In  Jjd  In  E)j^^ 

for  use  in  (3).  Thus,  given  a  field  model  from  which  to  calculate  (d  In  Byjd  In  E^),  we  can 

corrqiare  the  adiabatically  predicted  particle-flux  riKxlulation  with  the  riKxlulation  actually 

observed.  Alternatively,  we  can  use  the  CRRES  data  to  cortqwte  (e.g.,  as  a  function  of  L) 
the  value  of  (d  In  E^/d  In  E^)  that  would  be  required  in  order  to  account  adiabatically  for 

the  observed  modulation  of  with  E^.  In  the  present  work  we  analyze  (ZRRES  data  by 
each  of  these  methods. 


PRELIMINARY  RESULTS 

The  open  and  filled  circles  in  Figure  7  represent  stormtime  observations  of  the 
ratio  (A  In  /j^/A  In  B^)^^  taken  from  Orbit  306  (executed  1225-2219  UT,  28  November 

1990,  which  was  Day  332  of  the  year).  The  open  (filled)  circles  correspond  to  data  from 

the  outixrund  (inbound)  segment  of  Orbit  306.  The  meaning  of  the  numerator  in  this  ratio, 
which  constimtes  the  left-hand  side  of  (3),  is  the  following:  A  In  is  the  logaritiim  of  the 

ratio  of  the  values  of  observed  at  a  given  L  value  on  Orbit  306  and  inferred  from  the 
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best  fit  to  (cf.  Hgure  4,  solid  curve)  at  same  L  value  on  Orbit  297.  The  meaning  of  the 
denominator  is  as  follows:  A  In  is  the  logaridim  of  the  ratio  of  the  values  of  B,  as 

observed  at  CRRES  on  Orbit  306  and  as  given  by  the  Olson-Pfitzcr  model  along  Orbit  306 

(the  quiet-time  observations  of  B  on  Orbit  297  having  been  compatible  with  the  Olson- 
Pfitzer  model  at  that  time).  This  construction  is  intended  to  assure  that  A  In  0  only  as 

a  consequence  of  storm-associated  temporal  variations  in  the  magnetospheric  B  field  and 
not  (for  example)  as  a  consequence  of  magnetic-latitude  differences  between  different 
orbital  passes  of  the  spacecraft 

The  dashed  and  solid  curves  in  Bgure  7  represent  hypothetically  anticipated  values 
of  (A  In  JJHi  In  B^)jf^,  as  if  the  authentic  value  of  P  s  (d  In  B^/d  In  B^)  were  1.0  or  2.5, 

respectively.  As  we  have  noted  above,  the  case  P  =  1.0  would  correspond  to  a 
hypothetically  in  situ  variation  in  particle  energy  as  a  direct  response  to  the  local  change  in 
B,  whereas  the  case  P  -  2.5  would  correspond  (see  Appendix)  to  the  nonlocal  flux 
modulation  produced  by  temporal  variation  of  a  weak  but  patially  uniform  field 
perturbation  AB,  proportional  to  D^i  and  either  parallel  (or  antiparaUel)  to  the  geomag¬ 
netic  dipole  moment  Figure  7  shows  that  neidier  of  these  oversimplifications  comes  very 
close  to  describing  the  actual  situation. 

Thus,  we  invert  the  problem  and  ask  instead  how  p  s  (d  In  B^/d  In  B^)  would  have 
to  vary  with  B  in  order  to  reproduce  the  observed  (A  In  J^/A  In  For  this  purpose 

we  regard  (3)  as  a  linear  equation  to  be  solved  for  (d  In  B^/d  te  B^).  The  results  of  this 

operation  (applied  to  stormtime  data  from  the  outbound  segment  of  Orbit  306)  are  shown 

as  open  circles  (O)  in  Figure  8a.  As  thus  derived  from  the  CRRES  data,  the  required 
variation  of  the  parameter  (d  In  B^/d  In  B^)  with  R  (radial  distance,  measured  in  Bg)  is 

quite  systematic  and  easy  to  interpret  (at  least  qualitatively)  in  the  context  of  standard 
ring-current  models  [e.g.,  Hoffman  and  Bracken,  1967].  Results  of  the  same  operation, 
but  applied  to  stormtime  data  from  the  inbound  segment  of  Orbit  306,  are  shown  as  filled 
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drcles  (•)  in  Figure  8*.  This  variation  of  (d  In  In  with  R  is  more  difficult  to 
interpret  in  the  context  of  standard  ring-current  models. 

TENTATIVE  INTERPRETATION 

Our  Figure  9  constitutes  a  schematic  but  analytical  representation  of  the  equatorial 
magnetic-field  peiturbatioo  produced  by  the  stonntime  ring  current  It  is  based  largely  on 
Figures  2  and  3  of  Hoffman  and  Bracken  [1967],  showing  (the  component  of  AB 
perpendicular  to  the  equatorial  plane)  for  two  realizations  (corresponding  to  different 
values  of  of  their  optimal  “third-order”  model  Noteworthy  features  in  Figure  9  are 
the  inner  region  of  negative  at  values  of  R  ^  5.7  and  the  outer  region  of  positive  AB^ 

at  values  of  R  ^  5.7  (R  *  distance  from  dipole,  measured  in  Rg).  Since  the  positive 
fnagni»rir  flux  of  ling-cunent  origin  crossing  the  equator  at  R  ^  5.7  must  balance  the 
negative  magnetic  flux  of  ring-cuirent  origin  crossing  the  equator  at  R  ^  5.7,  the  effect  of 
the  ring  current  must  be  to  increase  the  value  of  R  at  which  any  magnetic  shell  (identified 
by  the  amount  of  magnetic  flux  enclosed)  crosses  the  magnetic  equator.  Such  “infiaiirai” 
of  a  magnetic  shell  [CoM/,  1966]  almost  invariably  reduces  the  value  of  Bq  (equatorial 
magnetic-field  intensity)  on  the  corresponding  field  lines,  even  if  the  equatorial  crossing 
points  now  occur  in  the  region  of  positive  AB^.  TTius,  the  value  of  (d  In  B^/d  In  B^)  should 

be  positive  (and  usually  >  1)  for  R  %  5.7,  negative  for  R  ^  5.7,  and  infinite  for  some  value 
of  R  »  5.7  (such  that  AB,  =  0).  This  is  essentiaUy  the  pattern  seen  in  Bgure  8n,  except 
that  the  singularity  in  P  appears  at  R  -  5.25  in  this  case. 

Equations  specifying  AB2  as  a  function  of  equatorial  R  for  the  ring-current  field 
model  shown  in  Hgure  9  are  given  in  the  Appendix.  We  regard  this  nrodel  as  protoQrpical 
of  the  actual  ring-current  field  and  not  merely  as  a  schematic  illustration.  The  radii  r,,  r2, 
rj,  r4,  and  (also  the  field  strengths  Bqi,  B23,  B4,  and  AB„„)  are  scalable  parameters  of 
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the  model  but  must  vary  in  unison  (maintaining  fixed  ratios)  in  order  to  conserve  magnetic 
flux.  Thus,  the  field  strengths  should  vary  in  proportion  to  and  the  radii  should  vary 
in  proportion  to  (being  reduced  by  a  factor  of  5.7/5.25  fiom  Rgure  9  to  Figure  8,  for 
example).  The  model  itself  is  defined  by  equation  (A6),  whereas  equations  (A7)-(A1 1) 
demonstrate  use  of  the  model  to  obtain  P  s  (d  In  B^Jd  In  B^)  as  a  function  of  equatorial  R. 

The  dashed  and  solid  curves  in  Hgure  8  are  derived  from  the  simplified  ring- 
current  model  described  in  the  Appendix.  The  solid  curve  corresponds  to  Bqi  = 
(2/3)£)st  (with  Djt  having  varied  from  about  -60  nT  to  about  -40  nT  during  Orbit  306). 
The  premise  here  is  that  the  ring  current  accounts  for  2/3  of  D^^,  while  currents 
consequently  induced  on  the  Earth’s  surface  for  the  other  third  of  [The  argument 
supporting  this  decomposition  is  quite  standard  in  geomagnetism:  The  model  ring-current 
immerses  the  Earth  in  a  nearly  uniform  magnetic  field  whose  intensity  we  call  Aqi  (as  it 
spans  a  region  extending  fiom  the  origin,  r  =  0,  to  a  radial  (fistance  r  =  rj  *  2  R^.  This 
field  would  be  excluded  from  actually  penetrating  the  Earth  if  it  induced  an  additional 
dipole  field  of  intensity  -Boi  at  the  poles  and  thus  of  intensity  at  the  equator.] 

The  curves  in  Rgure  8  correspond  to  as  if  there  were  no  such  induced 

currents.  Either  way,  the  model  described  in  the  Appendix  seems  to  account  remarkably 
well  for  the  pattern  of  data  points  in  Rgure  8a,  but  not  as  well  for  the  pattern  of  data 
points  in  Rgure  8b.  Thus,  we  regard  Rgure  8  with  mixed  emotions. 


PREDICTIVE  MODELING 

The  results  obtained  in  Figure  8a  suggest  that,  with  a  reliable  ring-current  field 
model  available,  it  might  be  possible  to  account  for  the  observed  stormtime  evolution  of 
radiation-belt  particle  intensities  directly  in  terms  of  variations  in  (i.e.,  without 
necessarily  invoking  data  from  an  onboard  spacecraft  magnetometer).  We  might  hesitate 
to  call  such  a  model  “predictive,”  since  there  is  usually  a  significant  delay  in  the  reporting 
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of  D^^  as  a  function  of  tiine.  However,  the  late  availability  of  Dj,  would  not  prevent 
energetic-electron  data  from  being  interpreted  eventually  in  the  context  of  a  tiine- 
dependent  ring  current.  In  any  case  the  electron  flux  should  respond  to  according  to 
the  relationship 

(A  In =  (1  -  [(y+ l)/2y](dlnJjJdlnE)f[j^}(dlnB^JdBQi) 

+  0  In  JJd  In  [id  In  B;/dBoi)  -  (d  In  B^/dBoj)] .  (6) 

This  result  is  obtained  from  (3)  upon  multiplication  of  both  sides  by  (d  In  B^ldB^O,  and  (6) 
is  implemented  by  identifying  Bq^  with  (2/3)Z>,|.  This  procedure  would  lead  to  an  explicit 
“prediction”  for  J^iE)  as  a  function  of  time  along  the  spacecraft  trajectory  if  means  were 
available  [e.g.,  Tsunitani  and  Baker,  1979]  for  predicting  the  future  ten^oral  variation  of 
Otherwise,  since  adiabatic  response  of  the  trapped  particle  population  to  changes  in 
ring-current  intensiiy  is  essentially  immediate,  we  would  still  have  in  (6)  a  means  of 
retrospectively  interpreting  stormtime  variations  in  energetic  particle  intensities  with 

Figure  10  illustrates  the  application  of  this  “predictive  imdeling”  capability  to  the 
670-keV  electron  data  fnnn  Orbit  306.  The  solid  curve  shows  the  quiet-time  radial  profile 
from  Orbit  297  (cf.  Hgure  4,  solid  curve).  The  dashed  curve  constitutes  an  adiabatic 
“prediction”  for  the  unidirectional  flux  along  the  outbound  segment  of  CRRES  Orbit 
306,  based  on  published  value  for  13—17  hr  UT  on  28  November  1990,  whereas  the 
open  (filled)  circles  represent  actual  outbound  (inbound)  observations  of  =  670 

keV  during  that  tune  interval.  For  the  purpose  of  predictive  modeling,  we  use  the 
representation  of  AB^  given  by  (A6)  in  the  Appendix  and  scale  this  by  setting  Bqi  = 
(2/3)Dsf  onboard  CRRES  magnetometer  for  this  purpose  except 

indirectly,  by  using  Figure  8a  to  set  =  5.25  (whereupon  r^,  r-^,  and  also  have 

values  7.9%  smaller  than  those  specified  in  the  Appendix). 
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It  may  be  possible  to  identify  a  systematic  relationship  between  the  {r„}  and  D^^ 
from  data  presently  available.  For  example,  more  negative  values  of  D^^  should 
correspond  to  deeper  penetrations  of  ring-current  ions  into  the  magnetosphere  and  thus  to 
smaller  values  for  the  {r„}.  For  this  reason  there  should  be  at  least  a  positive  correlation 
between  Dj,  and  the  {r„}  for  «  =  1,  2,  3,  4,  5.  Perfa^s  there  is  even  a  true  functional 
relationship.  If  so,  this  would  let  us  specify  the  {r„)  during  future  events  without  recourse 
to  in  situ  particle  and  field  data  in  each  case.  The  immediate  goal,  however,  is  to  learn 
whether  the  adiabatic  variation  of  a  realistic  ring-current  model  (even  with  some  of  its 
parameters  determined  by  the  local  data)  can  lead  in  practice  to  reliable  predictions  for  the 
stormtime  modulation  of  electron  flux  at  various  relativistic  energies. 

In  fact,  agreement  between  the  observed  and  predicted  stormtime  profiles  of  Jj^  in 
Hguie  10a  is  not  especially  good,  as  the  adiabatic  prediction  exceeds  the  observed 
stormtime  Jj^  by  0.4-1.4  orders  of  magnitude.  The  adiabatic  model  even  predicts  an 
increase  in  the  stormtime  electron  flux  (E  =  670  keV)  over  quiet-time  values  at  L  5.3  (R 
^  5.1  at  A.  >»  10°),  whereas  the  observations  (open  circles)  show  an  order-of-magnitude 
decrease  that  actually  deepens  somewhat  with  increasing  L.  These  initially  disappointing 
results  call  for  a  further  analysis  of  the  underlying  observational  data,  to  learn  what  went 
wrong  on  the  logical  path  from  Figure  8a  to  Figure  10a.  (Agreement  between  the  dashed 
curve  and  the  filled  circles  in  Figure  10b  is  surprisingly  good  by  comparison.) 

FURTHER  ANALYSIS 

A  partial  answer  to  the  question  of  what  went  wrong  between  Figure  8a  and 
Figure  10a  is  provided  in  Figure  11,  which  offers  separate  plots  of  A  In  7^  s^tl  A  In  R,  the 

observational  quantities  whose  ratio  appears  in  Figure  7.  Here  the  open  (filled)  circles 
represent  values  of  A  In  derived  from  data  acquired  on  the  outbound  (inbound)  segment 

of  Orbit  306,  whereas  the  +  (x)  signs  represent  values  of  A  In  R  derived  from  data 
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acquired  on  the  outbound  (inbound)  segment  of  the  same  orbit  Ihe  magnetic-&ld  data 
from  the  outbound  segment  of  Orbit  306  show  that  values  (+)  of  A  In  B  were  uniformly 

negative  (rather  than  negative  for  R  <  5.25  and  positive  foi  R  >  5.25).  This  pattern 
contradicts  the  premise  underlying  Figures  8  and  9,  whereby  the  value  of  (d  In  B^Jd  In  S^) 

is  supposed  to  have  become  singular  atr^r^  because  A  In  s  0  there.  It  seems  instead 
that  (d  In  B^^/d  In  B^)  became  singular  at  =  5.25  in  Hgure  8a  because  the  coefficient  of 
(d  In  Byjd  In  B^)  in  (3)  went  to  zero  there.  However,  the  coefficient  of  (d  In  B^Jd  In  B^)  in 

(3)  depends  entirely  on  the  quiet-time  radial  profile  and  energy  spectrum,  not  on  the 

structure  of  the  magnetic  storm.  Excellent  agreement  found  in  Rguie  8a  between  die 
observed  and  expected  profiles  of  (d  In  B^Jd  In  Bf^  thus  seems  in  retrospect  to  have  been 

largely  accidental  This  eiqietience  should  serve  as  a  warning  against  the  over- 
interpretation  of  observational  data. 

It  might  be  argued  that  die  outbound  segment  of  Orbit  306  was  executed  at  too 
high  a  magnetic  latitude  (10.7^  to  9.4**)  for  the  present  analysis  to  tqiply.  The  firfH  model 
described  in  the  Appendix  spedGes  AB  at  (and  thus  perpendicular  to)  the  equatorial  plane, 
whereas  AB  should  be  more  nearly  parallel  to  the  equatorial  jdane  at  latitudes  not  fru  from 
X  =  0**  [e.g.,  Hoffinan  and  Bracken,  1965].  If  magnetic  latitude  is  the  source  of  our 
troubles,  however,  then  die  need  to  account  for  nontrivial  conservation  of  the  second 
adiabatic  invariant  will  greatly  complicate  to  predict  the  electron-flux  modulatitms  from  a 
simple  model. 

We  should  not  have  expected  great  predictive  success  for  the  particle  and  field 
data  acquired  on  the  inbound  segment  of  Orbit  306.  Hguie  11  shows  that  A  te 

(represented  inbound  by  the  x  signs)  was  positive  where  it  should  have  been  negative 
(viz.,  at  ^  5.6)  and  negative  where  it  should  have  been  positive  (viz.,  atR^  5.7)  on  the 
inbound  segment,  althou^  was  very  well  behaved  throughout  (varying  from  -57  nT 
to  -50  nT  while  CRRES  was  inbound  on  Orbit  306).  The  anomalous  sign  of  A  In  in 
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Hguie  1 1  cannot  be  attributed  to  observations  made  at  too  high  a  magnetic  latitude, 

X  ranged  only  from  >6.2**  to  -1.0**  on  the  inbound  segment  of  Orbit  306.  This  was  as 
nearly  equatorial  a  pass  as  CRRES  ever  made  through  the  outer  radiation  belt 

Anyway,  Hguies  86  and  106  offer  an  analysis  of  Orbit  306  inbound  data  in  the 

format  of  Hguies  8a  and  10a.  The  filled  ciicles  in  Hguie  86  show  the  variation  that 
would  be  required  of  (d  In  In  in  order  to  account  adiabaticaUy  for  the  observed 

electron-flux  modulation.  The  solid  and  dashed  curves  show  the  radial  variation  that 
would  be  expected  of  (d  In  B^d  In  B^  if  Bqj  were  equal  to  D„  or  to  2/3  of 

respectively.  There  is  again  little  difference  between  the  dashed  and  solid  curves,  mainly 
because  the  field  of  the  induced  dipole  is  almost  negligible  at  B  4.  However,  there  is 

almost  no  resemblance  in  Hguie  86  between  these  curves  and  the  data  points  representing 
required  values  of  (d  In  B^/d  In  B^).  Hus  is  in  contrast  to  Hguie  8a,  in  which  the 

agreement  seemed  excellent  Hguie  106  shows  a  comparison  between  starmdme  values 
of  (jfiUed  dicles)  observed  inbound  on  Orbit  306  and  predictions  (dashed  curve)  for  the 

stormtime  based  on  (A6)-<A14),  as  iqrplied  to  the  quiet-imie  profile  (solid  curve  here 

and  in  Hguie  4)  and  to  fits  of  the  quiet-tinae  energy  spectra  observed  on  Orbit  297. 
Agreement  between  the  dashed  curve  and  the  data  points  for  while  still  not  excellent 

is  (as  noted  above)  actually  somewhat  better  in  Hguie  106  than  in  Hgure  10a.  Of  course, 
the  observed  profile  of  A  In  B^  in  Hguie  11  for  the  inbound  segment  of  Orbit  306 

disagrees  entirely  (as  noted,  even  in  sign)  from  what  Hguie  9  would  imply.  We  nurst 
conclude  that  adiabatic  response  to  a  simple  ring-current  fidd  nnodel  has  not  accounted 
well  for  the  energetic-electron  modulation  encountered  during  the  storm  of  27-28 
November  1990.  (We  lack  data  from  the  onboard  magnetometer  for  Orbit  30S,  during 
which  D^^  attained  -136  nT,  its  most  negative  value  of  this  storm.) 
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ANALYSIS  OF  A  LARGER  STORM 

The  results  described  above  came  from  a  moderate  magnetic  storm,  to  which  we 
had  expected  an  essendaOy  adiabatic  response  from  the  tipped  energetic  electrons.  The 
response  actually  observed  was  difficult  to  interpret,  at  least  partly  because  values  of  B 
(and  thus  of  A  In  R)  obtained  from  the  onboard  magnetometer  data  showed  quite  a 
different  variation  with  R  than  would  have  been  expected  from  standard  ring-current  field 
nuxiels.  (Perhaps  D^^  was  not  negative  enough  during  Orbit  306  to  provide  a  clean 
example  of  adiabatic  response.) 

Here  we  perform  a  similar  analysis  of  data  from  the  major  storm  of  9-10  July 
1991,  which  were  Days  190-191  of  the  year.  To  mtq)  flux  profiles  and  obtain  energy 
spectra  for  the  preceding  quiet  interval,  we  chose  Orbit  842  (executed  1449-0043  on  7-8 
July  1991).  The  spectra  encountered  on  Orbit  842  were  well  represented  as  exponential 
Their  ^-folding  energies  were  shown  in  Figure  6  above.  StormtitnB  data  used  in  oHiqnling 
Figure  12  (in  the  same  format,  except  for  ordinate  scales,  as  Hgure  1 1)  are  from  Orbit  846 
(executed  0625-1619  on  9  July  1991).  The  most  negative  value  (-198  nT)  of  attained 
during  this  storm  occurred  during  Hour  15  of  that  day.  As  Hgure  12  shows,  however, 
both  the  magnetic  field  and  the  670-keV  electron  flux  were  quite  unsettled  during  this  time 
interval,  to  the  extent  that  systematic  interpretation  ^rpears  difficult  in  this  case  also. 

We  have  nevertheless  tried  to  estimate  values  of  (d  In  R^/d  In  B^)  that  would  have 
been  required  in  order  to  account  adiabatically  for  observed  values  of  (d  In  JJb  In 

Tbf^  results  are  shown  in  Hgute  13  (same  format  as  m  Hgure  8).  Open  (filled)  circles 
correspond  to  outbound  finbound)  data  from  Orbit  846.  Predictions  for  (d  In  B^/d  In  B^), 

shown  as  solid  curves  for  Bq^  =  (2/3)Djt  and  as  dashed  curves  for  Bqi  =  differ  sii^tly 
between  the  inbound  and  outbound  passes.  These  minor  differences  arise  in  part  because 
of  hourly  changes  in  the  value  of  D^^,  which  appears  (via  Bqi)  in  the  numerator  and 
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denominator  of  (A14).  Nfoieover,  because  of  differences  in  ^Mcecraft  larinitfe  outbound 

and  inbound,  the  same  value  of  R  corresponds  to  different  values  of  L,  and  thus  to 
different  quiet-timc  values  of  £o  ^  ^  Id  between  the 

outbound  and  inbound  segments  of  Orbit  846.  Agreement  between  the  observed  and 
predicted  values  of  (d  In  B^d  In  is  not  good  on  the  outbound  pass  (Figure  13u!)  in  any 

case,  but  seems  quite  promising  early  on  the  inbound  pass  (Figure  13h). 

Hgure  14  (same  format  as  m  Figure  10)  shows  the  quiet-tune  flux  profile  (solid 
curves)  and  stormtime  predictions  (dashed  curves)  for  the  flux  profile  along  CRRES  Orbit 
846,  both  outbound  and  inbound.  The  predictions,  based  on  the  model  ling-cuirent  feM 

described  in  the  Appendix,  follow  from  (6)  for  figi  =  (2/3)D,i.  Open  (filled)  circles  in 
Figure  14  represent  actual  outbound  (inbound)  observations  of  on  Orbit  486.  There  is 

actually  reasonable  agreement  between  the  stormtime  predictions  and  the  data  points  at  R 
6,  as  well  as  with  some  of  the  isolated  data  points  at  lower  values  of  R,  Qmsiderable 
scatter  in  the  data,  however,  makes  it  difficult  to  be  more  precise  than  this,  and  many  of 
the  stormtime  data  points  still  fall  well  below  the  conespcMiding  stormtime  predictions. 

FURTHER  CONSIDERATIONS 

One  reason  why  we  are  consistently  overestimating  the  stormtime  electron  flux, 
even  by  an  order  of  magninide  in  several  instances,  is  that  our  basic  equations  require  us 

to  take  account  of  the  ‘Vadial  gradient”  of  the  quiet-time  electron  flux  profile  through  the 
factor  (d  In  JJd  In  in  (3)  and  (6).  Since  our  radial  profiles  tend  to  attain  tbeir 

maxima  at  L  »  4-5,  outward  displacement  of  the  profile  makes  a  positive  contribution  to 
A  In  /j,  over  most  of  the  range  of  R  shown  in  Figures  10  and  14.  This  tends  to  offset  the 
negative  contribution  to  A  In  Tj,  from  the  particle-energy  loss  associated  with  magnetic- 
shell  inflation  by  the  stoimtinie  ring  current 
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It  would  be  wrong,  of  course,  just  to  neglect  the  factor  (3  In  JJd  In  in  a 

serious  calculation.  However,  this  turns  out  to  be  an  instructive  exercise  nevertheless,  as 
it  shows  that  adiabatic  energy  loss  alone  is  usually  inadequate  to  account  for  the  typical 
stormtime  drop-out  of  energetic-electron  fluxes  at  L  5  even  if  the  accon^ying 
displacement  of  the  radial  profile  is  neglected.  In  order  to  sirrq)Iify  the  estimate  that 
supports  this  inference,  we  tty  it  at  r  =  r5  in  the  field  nnodel  defined  by  (A6).  This  has  the 
effect  of  making  =  0  at  the  location  of  interest,  which  roughly  corresponds  (since 
in  =  5.1a  in  Figure  9)  to  the  mean  r  value  at  which  the  electron-flux  observations  shown  in 

Figure  1  were  made.  Since  the  induced-dipole  term  specified  by  (A8)  is  negligible  at  such 
a  large  geocentric  distance,  we  can  neglect  the  term  (d  In  B^/cIBqj)  in  (A  12a),  substitute 

for  <l>re  there  by  invoking  (A7a),  and  finally  solve  (A12a)  as  a  differential  equation  to 

obtain 

AlnB^  »  —  31njl— 0.48(r5^/jiE)Boij  “  -3In|l-0.32(r5^/p,E)At|  (^) 

at  r  =  r^.  Our  quiet-time  estimates  for  (3  In  JJb  In  at  £  =  670  keV  in  this  region, 

obtained  from  exponential  fits  such  as  those  in  Figure  5,  are  listed  in  Table  3  for  the  15 
events  thus  analyzed.  From  these  fits  we  had  obtained  Eq  150—270  keV  at  L  *  5.5— 6.0, 
and  so  the  spectral  parameter  (3  In  JJd  In  £)jrx  ™  amounted  to  —3.5  ±  1.0  at 

E  =  670  keV  for  the  numbered  events  in  Table  1. 

Adiabatic  energy  loss  might  thus  reduce  the  trapped  electron  flux  by  a  factor  ~  2 
for  Djt  =  -35  nT  or  by  a  factor  ~  10-15  for  =  -135  nT.  (More  precise  numerical 
results  corresponding  to  the  15  analyzed  events  are  listed  in  Table  3.)  Reduction  factors 
actually  observed  range  from  about  10—50  for  Djj  =  -35  nT  to  about  100—400  for  D^^  =  — 
135  nT  (cf.  Table  1).  We  infer  that  adiabatic  energy  loss  alone  is  typically  inadequate  to 
account  for  the  entire  stormtime  drop-out  of  energetic-electron  fluxes  at  L  ^  5,  even  if  the 
offsetting  outward  motion  of  the  radial  profile  is  neglected.  For  several  events  in  Table  3, 
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however,  adiabatic  energy  loss  would  account  for  at  least  the  square  root  of  the  factor  by 
which  the  observed  quiet-time  flux  exceeded  the  observed  stormtimc  flux  at  £  =  670  kcV, 
and  thus  for  at  least  half  the  logarithmic  decrease.  Nevertheless,  we  must  conclude  that 
the  typical  stormtime  drop-out  in  lelativistic-etectron  flux  corresponds  at  least  in  pan  to  a 
genuine  loss  of  trapped  particles  from  the  magnetosphere. 

Such  an  inference  is  not  entirely  surprising,  as  relativistic  electron  precipitation 
(REP)  has  long  been  known  [Foitush  et  ai.,  1961]  and  understood  [Thorne  and  Kennel, 
1971]  to  accompany  the  main  phases  of  geomagnetic  storms.  The  present  work,  however, 
constitutes  a  careful  and  quantitative  study  of  the  relative  contributions  of  adiabatic 
response  and  electron  precipitation  to  the  observed  stormtime  modulation  of  relativistic 
electron  intensities.  In  this  context  it  does  seem  somewhat  surprising  that  even  tne  smaller 
storms  show  evidence  suggesting  a  significant  loss  of  relativistic  electrons  from  the 
magnetosphere.  Indeed,  three  of  the  five  largest  flux-reduction  factors  100)  in  the  last 
column  of  Table  3  correspond  to  storms  in  which  did  not  even  reach  -SO  nT. 

DISCUSSION  AND  SUMMARY 

Physical  processes  treated  in  the  present  work  are  essentially  adiabatic  and  thus 
essentially  kinematical  [cf.  Schulz  and  Lanzerotti,  1974,  p.  46].  However,  as  thei^  ^)pear 
to  account  only  partially  for  the  major  variations  observed  in  outer-zone  electron 
intensities  during  geomagnetic  storms,  we  can  infer  from  such  an  arfiahatir  study  that  a 
significant  part  of  the  stormtime  reduction  in  outer-zone  radiation  intensiQr  reflects  a 
genuine  loss  of  trapped  energetic  electrons  from  the  magnetosphere.  Thus,  even  a  purely 
adiabatic  study  such  as  this  can  lead  to  ureful  estimates  for  the  level  at  which  non- 
adiabatic  processes  are  operating  concurrently.  Standard  time-averaged  models  of  the 
radiation  environment  obviously  omit  stormtime  variability  altogether  and  so  inherently  fail 
to  distinguish  between  adiabatic  and  non-adiabatic  variations. 
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Dynamical  processes  involving  violation  of  one  or  more  of  the  adiabatic  invariants 
of  charged-particle  motion  are  presumed  to  accompany  the  adiabatic  processes,  especially 
under  stormtime  conditions.  Some  of  these  dynamical  processes  (e.g.,  wave-particle 
interactions)  lead  to  the  loss  of  trapped  particles,  while  others  (e.g.,  unsteady 
magnetospheric  convection)  allow  plasma-sheet  ions  and  electrons  to  reach  the  inner 
magnetosphere  and  thereby  establish  the  stormtime  ring  current.  Unsteady 
magnetospheric  convection  must  contribute  simultaneously  to  the  diffusive  radial  transport 
of  more  energetic  particles  (such  as  outn-belt  electrons)  among  drift  shells.  However, 
diffusive  processes  are  not  treated  expliddy  in  the  present  work,  and  ten^oral  variations 
of  the  axisymmetric  part  of  the  ring  current  (as  measured  by  the  index)  are  unlikely  to 
produce  much  diffusive  transport  of  radiation-belt  electrons  [cf.  Falthantmar,  1965]. 

While  radial  diffusion  must  play  an  inqxvtant  role  in  the  long-term  (e.g.,  week-to- 
week)  evolution  of  radiation-belt  particle  intensities,  results  obtained  in  the  present  work 
demonstrate  instead  the  means  for  constructing  a  time-dqiendent  model  of  the  trapped- 
electton  environment,  mduding  the  short-term  (e.g.,  hour-to-hour)  adiabatic  modulation 
associated  with  a  time-varying  D^.  Such  a  model  would  be  qrplicable  to  any  and  all  levels 
of  magnemspheric  distortion  (as  measured  by  the  index),  although  our  present 
investigation  is  limited  to  stornis  of  iiKxierate  intensity  so  as  to  exclude  catastrophic  events 
in  which  (for  example)  the  particles  of  interest  might  escape  from  the  magnetosphere 
along  suddenly  opened  field  lines  [e.g.,  Lyons  et  al.,  1989]. 

If  the  observed  modulation  of  energetic-electron  intensities  by  were  purely 
adiabatic,  this  would  appear  to  have  offered  a  parametric  neans  of  remotely  sensing  the 
global  distribution  of  the  Earth’s  ring  current,  even  on  time-scales  short  compared  to  the 
orbital  period  of  the  observing  spacecraft  In  other  wraris,  electron-flux  nnodulation 
anywhere  in  the  equatorial  magnetosphere  would  non-locally  reflect  the  adiabatic  influence 
of  ring-current  intensity  variations  elsewhere.  So  far,  however,  the  stormtime  modulation 
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of  nclativistic-electron  fluxes  appears  to  be  at  best  partially  adiabatic,  and  this  would 
certainly  complicate  the  use  of  energetic-electron  data  to  determine  the  parameters  (cf. 
Figure  9)  of  a  global  ring-current  model. 

Our  14-month  CRRES  data  base  initially  included  15  events  deemed  most 
amenable  to  analyses  of  the  sort  described  above.  These  are  the  cases  assigned  “event 
numbers”  in  Tables  1  and  3.  We  have  in  fact  investigated  several  additional  geomagnetic 

storms  from  October  1990  through  October  1991  (see  Table  1)  with  similarly  mixed 
results.  Sometimes  the  “required”  value  of  (d  in  fi^/d  In  B^)  deduced  fiom  the  data  agrees 

with  expectations  (as  in  Figure  8a),  sometimes  not  (as  in  Hgure  8b).  In  most  cases  of 
“agreement,”  however,  the  onboard  CJUiES  magnetometer  reports  a  AB(r)  that  is  quite 
different  from  the  expected  form  (cf.  Figure  9).  Sometimes  the  predicted  stormdme  flux 
profile,  obtained  by  adiabatically  transforming  quiet-time  spectra  and  radial  profiles 
through  use  of  the  model  for  AB(r)  shown  in  Bgure  9,  agrees  fairly  well  (as  in  Figure  10b) 
with  the  stormtime  data  points,  sometimes  not  (as  in  Hgure  10a).  Sometimes  the 
adiabatic  decrease  in  electron  energies  associated  with  the  stormtime  decrease  in 
(toward  more  negative  values)  accounts  for  half  or  more  of  the  observed  logarithmic 
decrease  in  relativistic  electron  flux  (as  for  Events  #1  and  #7  in  Table  3),  sometimes  not 
(as  for  the  other  events  in  Table  3). 

In  summary,  the  methods  developed  here  for  analyzing  CRRES  data  on  the 
stormtime  variation  of  relativistic-electron  intensities  yields  tantalizingly  promising  results 
on  occasion,  but  not  consistently  so.  The  reasons  for  this  inconsistency  are  unclear.  One 
possibility  suggested  by  the  present  analysis  is  that  the  magnetic  field  perturbation  AB(r) 
produced  by  the  stormtime  ring  current  has  a  more  complicated  spatial  structure  than  that 
envisioned  by  standard  ring-current  models.  Another  possibility  is  that  the  present 
methods  of  analysis,  being  based  on  the  premise  of  equatorially  mirroring  particles,  do  not 
adequately  account  for  the  adiabatic  bounce  motion  of  the  relativistic  electrons  of  interest. 
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Finally,  a  third  (and  rather  likely)  possibility  is  that  the  typical  stormtinie  drop-out  in 
relativistic-electron  flux  is  not  entirely  an  adiabatic  nxxlulation  but  corresponds  at  least  in 
part  to  a  genuine  loss  of  trapped  particles  from  the  magnetosphere.  If  so,  our  methods 
inherently  provide  a  quantitative  measurement  of  the  fraction  of  particles  thus  lost  in  each 
event  subsequently  analyzed. 
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APPENDIX:  REPRESENTATION  OF  RING-CURRENT  FIELD 

The  addition  of  a  unifoim  southward  (or  northward)  field  AB  =  zAfi^  parallel  (or 
anti-parallel)  to  the  geomagnetic  dipole  nK)nient  -li^z  would  lead  to  a  distortion  of  field 
lines  fi’om  the  usual  dipolar  configuration  specified  by  r = la  sin^fi,  where  a  (=  1  /?£)  is  the 
planetary  radius  and  L  is  a  dimensionless  label  inversely  proportional  [cf.  Roederer,  1970, 
p.  107]  to  the  magnetic  flux  <1>  =  enclosed  by  the  corresponding  magnefic  shell. 

The  equation  of  a  field  line  thus  distorted  [e.g..  Hill  and  Rassbach,  1975]  becomes 

r  =  Lfl[l-(r^Afl,/24E)]sin^e  ,  (Al) 

where  r  is  the  “geocentric”  distance  (actually  the  distance  from  the  point  dipole  to  the 
point  of  interest)  and  6  is  the  magnetic  colatitude.  Moreover,  the  value  of  L  in  (Al) 
remains  equal  to  The  magnetic  shell  of  interest  is  thus  “inflated”  (relative  to 

the  corresponding  dipolar  shell)  for  AB,  <  0  and  “cornpresscd”  for  AB^  >  0. 

The  field  model  that  leads  to  (Al)  is  not  a  good  representation  of  the  ring-current 
field.  Among  other  deficiencies,  it  is  current-free.  The  oversimplified  model  does, 
however,  conveniently  illustrate  the  issues  that  arise  in  nxxleling  ring-current  effects  on 
energetic  trapped-particle  distributions.  Most  notably  for  the  present  context,  the 
equatorial  magnetic  intensity  at  “geocentric”  distance  r  in  the  distorted  field  is  given  by 

Bo  =  =  (^E/'■^)[l  +  (^^AB,/HE)I.  (A2) 

whereas  the  equatorial  magnetic  intensity  on  the  dipolar  magnetic  shell  having  the  same  L 
value  (and  thus  enclosing  the  same  amount  of  magnetic  flux)  is 

Bo  =  (A3) 
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Thus,  a  particle  of  constant  third  invariant  O  sees  its  equatorial  B  value  (called  here) 

decrease  during  the  storm  by  a  factor  equal  to  the  ratio  between  the  right-hand  sides  of 
(A2)  and  (A3).  It  follows  that 


_  dln^  ^  (3/2)(rV4E) 

dAS,  “  dAB,  l-(l/2)(r^AB,/pE) 

and  that 

<‘1°^  _  (r^/He) 

dAB,  ■  l  +  (r’4fi,/|»E)  ■ 

whereupon 

_  (5/2)+(r[/iB,/iiE) 
dlnB;  "  l-(l/2)(rW,/fiB)  ' 


(A4a) 


(A4b) 


(A5) 


The  familiar  result  that  (d  In  B^/d  In  B^)  3  2.5  is  recovered  from  (A5)  by  taking  the  limit 
(/3AB,/He)-»0. 


To  illustrate  the  effect  of  a  mote  realistic  representation  of  the  ring-current  field  on 
geomagnetically  trapped  particles,  we  have  devised  an  analytical  representation  of  die 
fonn 


=  Bqi  ,  0  S  r  ^  ri 

(A6a) 

=  ^01  [(^“'■l)/(^2~n)K^23  “®0l)  »  '‘1  ^  ^  ^  ^2 

(A6b) 

~  ^23  *  ^2  ^  r  ^ 

(A6c) 

^rc 

-  B23  +  [(r-ri)Kr^-r^)](B^-B23),  ^  r  ^ 

(A6d) 

^rc 

=  ir^lr)\ir'^ -ri)l{rl -rl)]B^,  r  > 

(A6«) 
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for  the  north-south  (z )  of  the  equatorial  magnetic  field  produced  by  the  ring  current.  This 
representation  is  plotted  in  Figure  A  for  =  2.00,  =  3.59,  r-^  =  4.00,  =  4.9861 18279, 

and  =  5.70  (all  measured  in  Earth  radii,  /fg).  We  think  of  Bq^  (<  0)  as  2/3  of  D^^.  The 
other  third  comes  from  the  induced  dipole  (id),  which  diverts  ring-current  field  lines 
around  the  Earth.  To  the  other  “amplitude”  parameters  in  (A6)  we  assign  the  values  B23  = 
1.6i?oi  ^4  =  0.47228281 7^01  •  The  values  of  and  arc  determined  by  requiring 

that  (A6d)  and  (A6e)  join  smoothly  at  r  =  subject  to  the  constraint  that  the  maximum 
positive  AB^  (attained  at  r  =  be  equal  to  -0.16Boi-  The  other  parameters  in  (A6) 
are  assigned  on  the  basis  of  visual  impressions  of  various  ring-current  field  models  [e.g., 
Hoffman  and  Bracken,  1967]  and  adjusted  as  necessary  so  as  to  conserve  magnetic  flux. 
(The  positive  flux  at  r  >  r5  balances  the  negative  flux  at  r  <  r5  to  about  one  part  in  10^  in 
this  example.)  The  model  ring-current  field  in  Figure  A  scales  in  field  intensity  with  D^^ 
and  can  be  scaled  radially  by  applying  a  common  factor  to  all  the  {r„}  for  n  =  1,  2,  3, 4,  5. 
Simulations  of  the  particle  transport  needed  to  produce  the  stormtime  ring  current  [Chen 
et  al.,  1994]  suggest  that  values  assigned  to  the  {r„}  should  vary  inversely  with  |D,J. 

From  the  equatorial  field  model  specified  by  (A6),  it  is  easy  to  calculate  the 
amount  of  magnetic  flux  outside  any  magnetic  sheU  of  equatorial  radius  r.  The  results 
needed  for  use  in  the  present  study  are 

=  -1.287tr52(r5/r)[l-l(r5/r)^]Bbi  ,  r  ^  (A7a) 


and 

<I>re  =  (rt  /  3)[(B4  -  B23 )  /  ('4  -  '3  )1(2'4  -  ) 

-  1.287cr5^(r5/r)[l— i-(r5/r)^]floi  +  itir^ -r^)B23  ,  ^  r  ^  {Mb) 

In  particular,  the  amount  of  positive  flux  at  r  >  r5  is  equal  to  -0.%7tr5  Bqi- 
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The  magnetic  field  generated  by  the  ring  current  is  {approximately  uniform,  not 
only  at  the  equator  for  r  <  rj  *  2  R^,  but  throughout  a  spherical  volume  of  this  radius 
surrounding  the  Earth  [e.g.,  Hoffman  and  Bracken,  1965].  Under  the  {qjproximation  that 
this  magnetic  field  is  excluded  firom  penetrating  the  Earth  itself,  the  ring  current  can  be 
considered  to  “induce”  a  dipole  of  moment  (fl3/2)Boi.  which  is  generated  by  an  azimuthal 
current  (in  a  direction  opposite  to  the  ring  current)  on  the  Earth’s  surface.  The  (z) 
component  of  this  induced-dipole  (id)  field  perpendicular  to  the  equatorial  plane  is  given 
by 

ABid  =  -Boi.  0  <.  r  <  a;  ABjd  =  (.aV2r^)BQi,  r  >  a.  (A8) 

and  thus  contributes  1/3  of  the  equatorial  that  we  identify  with  The  other  2/3  of 
Djt  comes  from  (A6a). 

The  total  magnetic  flux  exterior  to  a  magnetic  shell  of  equatorial  radius  r  is  thus 
equal  to 

<I>  =  (7t/r)(2pE  +  a^Boi)  +  (A9) 

with  On;  given  by  (A7)  for  r  ^  ry  The  equatorial  magnetic  intensity  on  the  dipolar 
magnetic  shell  having  the  same  L  value  (and  thus  enclosing  the  same  amount  of  magnetic 
flux)  is 

Bo  =  +  +  (AlO) 

whereas  the  equatorial  magnetic  intensity  at  “geocentric”  distance  r  in  the  distorted  field  is 
equtil  to 

£o  =  =  (^E/^^Xl  +  (a^Bol/2pE)  +  (^^ABn./pE)] .  (AH) 
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with  given  by  (A6d)  or  by  (Abe),  whichever  is  aijpropriate,  for  r  ^  r^.  Thus,  a 
particle  of  constant  third  invariant  <1>  sees  its  equatorial  B  value  (called  here)  decrease 

during  the  storm  by  a  factor  equal  to  the  ratio  between  the  right-hand  sides  of  (All)  and 
(A  10).  It  follows  that 


_  (3aV24E)-K3r<b^/2icflh,4p) 
l  +  (fl  A)i/24e)  +  (^ic/2JI|1e) 


(A  12a) 


and  that 

_  (a^/ 2pe)  +  (r^A^/ 
dfloi  1  +  (^i/2Pe 

whereupon 


(A12b) 


<n°^  ,  ,  _  (3/2)«^„[H^(■^A^/^B)l  „  ,3, 

dtaflj  io^AB;,[l  +  (f«„/2iiHB)]  ■ 

since  the  induced-dipole  term  makes  tally  a  iKgligible  contribution  to  AO  at  r  ^  r3  and 
shifts  only  slightly  (from  r  »  r^)  the  value  of  r  at  which  AB^  =  0.  Thus,  for  example,  we 
should  expect  to  find 


dln^  _  _  3tl - (l/4)(«i/r)’l(l-0.64(r? /4B)flo,[l -Wr)^]) 

dlnB;  ”  [1  -  (rj/  r)’lU-0.64(r/  /HE)£b,[l  -  (1/  dXr,/  r)’]) 

(with  Pe  =  30500  nT-/?E)  5.  We  have  used  this  last  result,  as  specified  1^ 

(A13),  for  generating  the  dashed  and  solid  curves  in  Hgure  8. 
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HGURE  CAPTIONS 

Rgurc  1.  Omnidirectional  electron  fluxes  in  two  energy  channels,  averaged  over  the  time 
spent  by  CRRES  between  L  =  4.5  and  L  =  6.5  on  each  orbit  For  comparison:  Histogram 
showing  hourly  values  of  D^^  for  the  same  50  days  (27  October  through  6  December 
1990). 

Hgure  2.  Local  (+)  and  equatorial  (□)  pitch-angle  distributions  of  electrons  in  the  561- 
958  keV  energy  channel  during  two  partial  orbits  executed  28-29  November  (Days  332- 
333)  1990.  This  was  late  in  the  recovery  phase  of  a  storm  during  which  had  (in  Hour 
22  of  Day  331)  reached  -136  nT. 

Figure  3.  Local  (+)  and  equatorial  (□)  pitch-angle  distributions  of  electrons  in  the  561- 
958  keV  energy  channel  while  CRRES  was  at  L  >  6  on  Orbit  702,  executed  2031-0624 
UT  on  9-10  May  (Days  129-130)  1991.  CRRES  attained  apogee  during  the  second  hour 
of  Day  130.  This  was  during  the  recovery  fiom  a  small  storm  in  which  had  (at  Hour 
22  of  Day  129)  reached  -23  nT  as  its  extremum.  The  feature  of  particular  interest  here  is 
die  “butterfly”  form  of  the  pitch-angle  distribution. 

Figure  4.  Radial  profile  of  electron  flux,  as  observed  in  the  £  =  561-958  keV  energy 
rhannH  on  Orbit  297,  executed  1946-0537  UT,  24-25  November  1990  (Days  331-332 
of  the  year).  Dashed  (dotted)  curves  are  fits  to  outbound  (inbound)  data,  represented  by 
open  (filled)  circles,  respectively.  Solid  curve  is  a  joint  fit  to  both  outbound  and  inbound 
data.  Maximum  in  solid  curve  at  £  «  390  nT  corresponds  to  L  »  4.4,  as  CRRES  has  a 
low-inclination  orbit 

Figure  5.  Representative  quiet-time  energy  spectra  of  electrons  encountered  at  L  »  4.5 
and  L  =  5.5  on  Orbit  297,  both  outbound  and  inbound.  Symbol  Eq  denotes  characteristic 
(^-folding)  energy  of  exponential  fit  in  each  case. 


65 


43 


Hguic  6.  Characteristic  (^-folding)  energies  of  exponential  spectra  acquired  during  Orbit 
842  (executed  1449-0043  on  7-8  July  1991,  which  were  Days  188-189  of  the  year). 
Open  (filled)  circles  fitted  by  dashed  (solid)  parabola  correspond  to  outbound  (inbound) 
segment  of  Orbit  842. 

Rgurc  7.  Open  (filled)  circles  represent  left-hand  side  of  (3)  according  to  data  fi-om 

outbound  (inbound)  segment  of  Orbit  306.  For  comparison,  dashed  and  solid  curves 
represent  right-hand  side  of  (3)  for  |3  s  (d  In  B^/d  In  s  1.0  and  P  s  2.5,  respectively, 

with  energy  spectrum  and  radial  profile  based  on  data  acquired  during  Orbit  297. 

Figure  8.  Open  (filled)  circles  indicate  values  required  of  (d  In  Byjd  In  in  order  that 
corresponding  outbound  (inbound)  observations  of  ratio  (A  In  JJd  In  Bj^  on  Orbit  306 
satisfy  (3).  Solid  and  dashed  curves  represent  values  of  (d  In  B^ld  In  B^^  expected  fi'om 

(A13)  with  5oi  =  (2/3)Z),t  and  Bqi  respectively. 

Figure  9.  Ring-current  field  model  defined  by  (A6).  Radii  (n  =  1,  2,  3,  4,  5)  and  field 
strengths  (figi,  ^23*  ^  scalable  (but  not  independently  adjustable)  parameters 

of  the  model  (constrained  by  flux-conservation  requirement). 

Figure  10.  Open  (filled)  circles  represent  stormtime  electron  fluxes,  670  keV,  measured 
outbound  (inbound)  on  Orbit  306.  Dashed  curves  represent  stormtime  flux  predictions 
based  on  (6).  Solid  curve  represents  quiet-time  flux  profile  fi’om  Orbit  297. 

Figure  11.  Open  (filled)  circles  represent  stormtime  electron  fluxes,  670  keV,  measured 
outbound  (inbound)  on  Orbit  306,  relative  to  quiet-time  electron  fluxes  measured  on  Orbit 
297.  For  comparison:  Values  of  B  =  |B|  obtained  from  outbound  (+)  and  inbound  (x) 
measurements  on  Orbit  306,  relative  to  Olson-Pfitzcr  model  B  on  Orbit  306. 
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Hguie  12.  Open  (filled)  aides  represent  sttvmtiinB  election  fluxes,  670  keV,  measured 
outbound  (inbound)  on  Orbit  486,  relative  to  quiet-tune  election  fluxes  measured  on  Orbit 
482.  For  comparison:  Values  of  B  =  iB|  obtained  from  outbound  (4-)  and  inbound  (x) 
measurements  on  Orbit  486,  relative  to  Olson-Pfitzer  model  B  on  Orbit  486. 

Figure  13.  Open  (filled)  aides  indicate  values  required  of  (d  In  By/d  In  B(^)  in  order  that 
corresponding  outbound  (inbound)  observations  of  ratio  (A  In  JJd  In  B^)  on  Orbit  486 
satisfy  (3).  Solid  and  dashed  curves  renresent  values  of  (d  In  B^Jd  In  B^)  expected  from 

(A13)  with  Boi  =  and  ^oi  =  respectively. 

Figure  14.  Open  (filled)  aides  represent  stormtime  electron  fluxes,  670  keV,  measured 
outbound  (inbound)  on  Orbit  486.  Dashed  curves  represent  stormtime  flux  predictions 
based  on  (6).  Solid  curve  represents  quiet-time  flux  profile  from  Orbit  482. 
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Table  1.  Magnetic  Storms  Reinesented  in  Base  for  Present  Study 


Event 

No. 

Orbit 

Numbers 

1 

064-071 

2 

114-119 

3 

131-136 

4 

185-191 

218-222 

5 

232-240 

6 

257-265 

275-278 

7 

303-309 

318-322 

8 

340-344 

9 

464-472 

580-585 

586-600 

11 

663-668 

12 

701-707 

13 

710-716 

14 

720-725 

763-772 

15 

780-790 

825-828 

852-860 

900-903 

1042-1045 

Days  of  Year, 
1990  or  1991 


233.1- 236.4  (90) 

253.7- 256.1  (90) 
260.6-263.1  (90) 

282.8- 285.7  (90) 

296.4- 298.4  (90) 

302.1- 305.8  (90) 

312.4- 316.1  (90) 

319.8- 321.4  (90) 

331.3- 334.6  (90) 

337.4- 339.5  (90) 

346.5- 348.5  (90) 

356.4- 360.0  (90) 
079.9-082.4  (91) 
082.4-088.1  (91) 

113.9- 116.3(91) 

129.4- 132.3  (91) 

133.1- 136.0(91) 

137.2- 139.7  (91) 

154.8- 159.0  (91) 

161.8- 166.4  (91) 

181.1- 182.8(91) 

192.6- 196.5  (91) 

213.2- 214.9  (91) 

274.2- 275.9  (91) 


Most  Negative 
Value  of  Djt 


-98  nT 
-35  nT 
-47  nT 
-133  nT 
-65  nT 
-50  nT 
-28  nT 
-54  nT 
-136  nT 
-39  nT 
-41  nT 
-50  nT 
-40  nT 
-298  nT 
-31  nT 
-23  nT 
-70  nT 
-103  nT 
-219  nT 
-138  nT 
-56  nT 
-185  nT 
-111  nT 
-162  nT 


Corresponding 
Hour  &  Date 


3x,  Aug  22-23 
Hr  19,  Sept  10 
Hr  16,  Sept  18 
Hr  09,  Oct  10 
Hrs  5-6,  Oct  24 
Hr  03,  Oct  30 
3x,Nov  10-11 
Hr  12,  Nov  16 
Hr  22,  Nov  27 
Hr  10,  Dec  04 
Hr  06,  Dec  13 
Hr  23,  Dec  24 
Hr  23,  Mar  21 
Hr01,Mar25 
Hr  01,  Apr  25 
Hr  22,  May  09 
Hr  18,  May  14 
Hr  10,  May  17 
Hr  20,  June  05 
Hr  07,  June  11 
Hr  22,  June  30 
Hr  16,  July  13 
Hr  16,  Aug  02 
Hrs  3-4,  Oct  02 


Rux  Reduction 
Factor,  670  keV 


60 

50 

160 

380 

160 

12 

460 

90 

70 

310 

26 

60 

100 

20 

17 

6 

60 

80 

280 

250 

36 

100 

100 

280 


Note:  Most  negative  value  of  was  attained  three  times  each  during  Events  #1  and  #6. 
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Table  2.  SEP  Energy  Ranges  (and  Channel  Widths),  MeV 


Mode 

Sensor  A 

Sensor  B 

Sensor  C 

ELECl 

0.042-0.324  (0.0235) 

0.042-0.336  (0.0245) 

0.041-0.313  (0.0227) 

ELEC2 

0.164-4.93 

(0.397) 

0.171-5.12 

(0.413) 

0.170-5.11 

(0.412) 

PROTl 

0.875-6.60 

(0.478) 

0.916-6.70 

(0.482) 

0.920-6.80 

(0.490) 

PROT2 

2.5-38.7 

(3.01) 

2.2-33.7 

(2.62) 

2.0-30.4 

(2.37) 

PROT3 

35.8-80.2 

(3.70) 

31.2-69.9 

(3.22) 

28.2-63.1 

(2.91) 

PROT4 

45.0-94.0 

(4.08) 

45.0-105.0 

(5.00) 

45.0-110.0 

(5.42) 

ALPHA 

6.8-24.0 

(1.43) 

6.9-24.3 

(1.45) 

7.0-24.6 

(1.47) 
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Table  3.  Flux-Reduction  Factors:  Expected  (from  Adiabatic  Energy  Loss)  vs  Observed 


Event 

No. 

Orbit 

Numbers 

Spectral  Index, 

(0  In  JjJd  In  £)jr^ 

Most  Negative 
Value  of 

Expectd  Reductn 
Factor,  670  keV 

ObsTvd  Reductn 
Factor,  670  keV 

1 

064-071 

-4.40 

-98  nT 

8.8 

60 

2 

114-119 

-3.80 

-35  nT 

2.1 

50 

3 

131-136 

-3.63 

-47  nT 

2.6 

160 

4 

185-191 

-3.75 

-133  nT 

12.7 

380 

5 

232-240 

-3.20 

-50  nT 

2.5 

12 

6 

257-265 

-3.80 

-28  nT 

1.8 

460 

7 

303-309 

-3.10 

-136  nT 

9.6 

70 

8 

340-344 

-3.30 

-^1  nT 

2.2 

26 

9 

464^72 

-5.20 

-50  nT 

7.0 

60 

10 

580-585 

-4.30 

-40  nT 

2.5 

100 

11 

663-668 

-3.35 

-31  nT 

1.8 

17 

12 

701-707 

-3.60 

-23  nT 

1.3 

6 

13 

710-716 

-3.30 

-70  nT 

3.6 

60 

14 

720-725 

-2.60 

-103  nT 

2.0 

80 

:5 

780-790 

-3.24 

-138  nT 

10.7 

250 
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